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ABSTRACT 



A review is given on the consistency checks of Grand Unified Theories (GUT), which unify the 
electroweak and strong nuclear forces into a single theory. Such theories predict a new kind of 
force, which could provide answers to several open questions in cosmology. The possible role of 
such a "primeval" force will be discussed in the framework of the Big Bang Theory. 
Although such a force cannot be observed directly, there are several predictions of GUT's, 
which can be verified at low energies. The Minimal Super symmetric Standard Model (MSSM) 
distinguishes itself from other GUT's by a successful prediction of many unrelated phenomena 
with a minimum number of parameters. 

Among them: a) Unification of the couplings constants; b) Unification of the masses; c) Existence 
of dark matter; d) Proton decay; e) Electroweak symmetry breaking at a scale far below the 
unification scale. 

A fit of the free parameters in the MSSM to these low energy constraints predicts the masses of 
the as yet unobserved superpartners of the SM particles, constrains the unknown top mass to 
a range between 140 and 200 GeV, and requires the second order QCD coupling constant to be 
between 0.108 and 0.132. 



(Published in "Progress in Particle and Nuclear Physics, 33 (1994) 201".) 



^Email: Wim.de. Boer@cern.ch 
Based on lectures at the Herbstschule Maria Laach, Maria Laach (1992) and the Heisenberg- 
Landau Summerschool, Dubna (1992). 



"The possibility that the universe was generated from noth- 
ing is very interesting and should be further studied. A most 
perplexing question relating to the singularity is this: what 
preceded the genesis of the universe? This question appears 
to be absolutely methaphysical, but our experience with meta- 
physics tells us that metaphysical questions are sometimes 
given answers by physics. " 
A. Linde (1982) 
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Chapter 1 
Introduction 



The questions concerning the origin of our universe have long been thought of as 
metaphysical and hence outside the realm of physics. 

However, tremendous advances in experimental techniques to study both the 
very large scale structures of the universe with space telescopes as well as the 
tiniest building blocks of matter - the quarks and leptons - with large accelera- 
tors, allow us "to put things together" , so that the creation of our universe now 
has become an area of active research in physics. 

The two corner stones in this field are: 

• Cosmology, i.e. the study of the large scale structure and the evolution 
of the universe. Today the central questions are being explored in the 
framework of the Big Bang Theory (BBT)[|| |^, H, ||, |^, which provides a 
satisfactory explanation for the three basic observations about our universe: 
the Hubble expansion, the 2.7 K microwave background radiation, and the 
density of elements (74% hydrogen, 24% helium and the rest for the heavy 
elements) . 

• Elementary Particle Physics, i.e. the study of the building blocks of 
matter and the interactions between them. As far as we know, the build- 
ing blocks of matter are pointlike particles, the quarks and leptons, which 
can be grouped according to certain symmetry principles; their interactions 
have been codified in the so-called Standard Model (SM)0. In this model 
all forces are described by gauge field theories^, which form a marvelous 
synthesis of Symmetry Principles and Quantum Field Theories. The latter 
combine the classical field theories with the principles of Quantum Mechan- 
ics and Einstein's Theory of Relativity. 

The basic observations, both in the Microcosm as well as in the Macrocosm, are 
well described by both models. Nevertheless, many questions remain unanswered. 
Among them: 

• What is the origin of mass? 

• What is the origin of matter? 

• What is the origin of the Matter-Antimatter Asymmetry in our universe? 
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Figure 1.1: The evolution of the universe and the energy scale of some typical 
events: above the Planck scale of 10^^ GeV gravity becomes so strong, that 
one cannot neglect gravity implying the need for a "Theory Of Everything" to 
describe all forces. Below that energy the well known strong and electroweak 
forces are assumed to be equally strong, implying the possibility of a Grand Uni- 
fied Theory (GUT) with only a single coupling constant at the unification scale. 
After spontaneous symmetry breaking the gauge bosons of this unified force be- 
come heavy and "freeze out" . The remaining forces correspond to the well known 
SU {3)c ® SU (2) L (g) U{1)y symmetry at lower energies with their coupling con- 
stants changing from the unified value at the GUT scale to the low energy values; 
this evolution is attributed to calculable radiative corrections. Future accelera- 
tors are expected to reach about 15 TeV corresponding to a temperature of 10^^ 
K, which was reached about 10~^^ s after the "Bang". At about 10^ GeV the 
gauge bosons of the electroweak theory "freeze out" after getting mass through 
spontaneous symmetry breaking and only the strong and electromagnetic force 
play a role. About three minutes later the temperature has dropped below the 
nuclear binding energy and the strong force binds the quarks into nuclei (nucle- 
osynthesis). Most of the particles annihilate with their antiparticles into a large 
number of photons after the photon energies become too low to create new par- 
ticles again. After about hundred thousand years the temperature is below the 
electromagnetic binding energies of atoms, so the few remaining electrons and 
protons, which did not annihilate, form the neutral atoms. Then the universe 
becomes transparent for electromagnetic radiation and the many photons stream 
away into the universe. The photons released at that time are now observed as 
the 3 K microwave background radiation. Then the neutral atoms start to cluster 
slowly intostars and galaxies under the influence of gravity. 



• Why is our universe so smooth and isotropic on a large scale? 

• Why is the ratio of photons to baryons in the universe so extremely large, 
on the order of 10^''? 

• What is the origin of dark matter, which seems to provide the majority of 
mass in our universe? 

• Why are the strong forces so strong and the electroweak forces so weak? 

Grand Unified Theories (GUT)[^ |^, in which the known electromagnetic, weak, 
and strong nuclear forces are combined into a single theory, hold the promise of 
answering at least partially the questions raised above. For example, they explain 
the different strengths of the known forces by radiative corrections. At high en- 
ergies all forces are equally strong. The Spontaneous Symmetry Breaking (SSB) 
of a single unified force into the electroweak and strong forces occurs in such the- 
ories through scalar fields, which "lock" their phases over macroscopic distances 
below the transition temperature. A classical analogy is the build-up of the mag- 
netization in a ferromagnet below the Curie-temperature: above the transition 
temperature the phases of the magnetic dipoles are randomly distributed and 
the magnetization is zero, but below the transition temperature the phases are 
locked and the groundstate develops a nonzero magnetization. Translated in the 
jargon of particle physicists: the groundstate is called the vacuum and the scalar 
fields develop a nonzero "vacuum expectation value". Such a phase transition 
might have released an enormous amount of energy, which would cause a rapid 
expansion ("inflation") of the universe, thus explaining simultaneously the origin 
of matter, its isotropic distribution and the flatness of our universe. 

Given the importance of the questions at stake, GUT's have been under in- 
tense investigation during the last years. 

The two directly testable predictions of the simplest GUT, namely 

• the finite lifetime of the proton 

• and the unification of the three coupling constants of the electroweak and 
strong forces at high energies 

turned out to be a disaster for GUT's. The proton was found to be much more 
stable than predicted and from the precisely measured coupling constants at the 
new electron-positron collider LEP at the European Laboratory for Elementary 
Particle Physics -CERN- in Geneva one had to conclude that the couplings did 



not unify, if extrapolated to high energies [|10|, \12^ . 

However, it was shown later, that by introducing a hitherto unobserved sym- 
metry, called Supersymmetry (SUSY)|jT3|, |l4l, into the Standard Model, both 



problems disappeared: unification was obtained and the prediction of the proton 
life time could be pushed above the present experimental lower limit! 

The price to be paid for the introduction of SUSY is a doubling of the number 
of elementary particles, since it presupposes a symmetry between fermions and 
bosons, i.e. each particle with even (odd) spin has a partner with odd (even) 
spin. These supersymmetric partners have not been observed in nature, so the 
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only way to save Supersymmetry is to assume that the predicted particles are 
too heavy to be produced by present accelerators. However, there are strong 
theoretical grounds to believe that they can not be extremely heavy and in the 
minimal SUSY model, the lightest so-called Higgs particle will be relatively light, 
which implies that it might even be detectable by upgrading the present LEP 
accelerator. But SUSY particles, if they exist, should be observable in the 
next generation of accelerators, since mass estimates from the unification of the 



precisely measured coupling constants are in the TeV region |TT[ and the lightest 
Higgs particle is expected to be of the order of M^, as will be discussed in the 
last chapter. 

It is the purpose of the present paper to discuss the experimental tests of 
GUT's. The following experimental constraints have been considered: 

• Unification of the gauge coupling constants; 

• Unification of the Yukawa couplings; 

• Limits on proton decay; 

• Electroweak breaking scale; 

• Radiative 6 — > 57 decays; 

• Relic abundance of dark matter. 

It is surprising that one can find solutions within the minimal SUSY model, which 
can describe all these independent results simultaneously. The constraints on the 
couplings, the unknown top-quark mass and the masses of the predicted SUSY 
particles will be discussed in detail. 

The paper has been organized as follows: In chapters ^ to |^ the Standard 
Model, Grand Unified Theories (GUT) and Supersymmetry are introduced. In 
chapter ^ the problems in cosmology will be discussed and why cosmology "cries" 
for Supersymmetry. Finally, in chapter ^ the consistency checks of GUT's through 
comparison with data are performed and in chapter ^ the results are summarized. 
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Chapter 2 

The Standard Model. 



2.1 Introduction. 

The field of elementary particles has developed very rapidly during the last two 
decades, after the success of QED as a gauge field theory of the electromagnetic 
force could be extended to the weak- and strong forces. The success largely 
started with the November Revolution in 1974, when the charmed quark was 
discovered simultaneously at SLAC and Brookhaven, for which B. Richter and 
S.S.C Ting were awarded the Nobel prize in 1976. This discovery left little doubt 
that the pointlike constituents inside the proton and other hadrons are real, 
existing quarks and not some mathematical objects to classify the hadrons, as 
they were originally proposed by Gellman and independently by ZweigQ. 

The existence of the charmed quark paved the way for a symmetry between 
quarks and leptons, since with charm one now had four quarks {u, d, c and s) 
and four leptons (e, /i, i^e and u^), which fitted nicely into the SU{2) ® U{1) 
unified theory of the electroweak interactions proposed by Glashow, Salam and 
Weinberg (GSW) [H for the leptonic sector and extended to include quarks as 
well as leptons by Glashow, Iliopoulis and Maiani (GIM) [|1^ as early as 1970. 
Actually, from the absence of flavour changing neutral currents, they predicted 
the charm quark with a mass around 1-3 GeV and indeed the charmed quark 
was found four years later with a mass of about 1.5 GeV. This discovery became 
known as the November Revolution, mentioned above. 

The unification of the electromagnetic and weak interactions had already 
been forwarded by Schwinger and Glashow in the sixties. Weinberg and Salam 
solved the problem of the heavy gauge boson masses, required in order to explain 
the short range of the weak interactions, by introducing spontaneous symmetry 
breaking via the Higgs-mechanism. This introduced gauge boson masses without 
explicitly breaking the gauge symmetry. 

The Glashow- Weinberg-Salam theory led to three important predictions: 

• neutral currents, i.e. weak interactions without changing the electric charge. 
In contrast to the charged currents the neutral currents could occur with 

^Zweig called the constituents "aces" and believed they really existed inside the hadrons. 
This belief was not shared by the referee of Physical Review, so his paper was rejected and 
circulated only as a CERN preprint, albeit well-known p5[|. 
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leptons from different generations in the initial state, e.g. i/^e v^e 
through the exchange of a new neutral gauge boson. 

• the prediction of the heavy gauge boson masses around 90 GeV. 

• a scalar neutral particle, the Higgs boson. 

The first prediction was confirmed in 1973 by the observation of scattering 
without muon in the final state in the Gargamelle bubble chamber at CERN. 
Furthermore, the predicted parity violation for the neutral currents was observed 
in polarized electron-deuteron scattering and in optical effects in atoms. These 
successful experimental verifications [0] led to the award of the Nobel prize in 1979 
to Glashow, Salam and Weinberg. In 1983 the second prediction was confirmed 
by the discovery of the W and Z bosons at CERN in pp collisions, for which C. 
Rubbia and S. van der Meer were awarded the Nobel prize in 1985. 

The last prediction has not been confirmed: the Higgs boson is still at large 
despite intensive searches. It might just be too heavy to be produced with the 
present accelerators. No predictions for its mass exist within the Standard Model. 
In the supersymmetric extension of the SM the mass is predicted to be on the 
order of 100 GeV, which might be in reach after an upgrading of LEP to 210 
GeV. These predictions will be discussed in detail in the last chapter, where a 
comparison with available data will be made. 

In between the gauge theory of the strong interactions, as proposed by Fritzsch 
and Gell-Mann ITT], had established itself firmly after the discovery of its gauge 



field, the gluon, in 3-jet production in e e~ annihilation at the DESY labora- 
tory in Hamburg. The colour charge of these gluons, which causes the gluon 
self-interaction, has been established firmly at CERN's Large Electron Positron 
storage ring, called LEP. This gluon self-interaction leads to asymptotic freedom, 
as shown by Gross and Wilcek |1^ and independently by Politzer [0, thus ex- 
plaining why the quarks can be observed as almost free pointlike particles inside 
hadrons, and why they are not observed as free particles, i.e. they are confined 
inside these hadrons. This simultaneously explained the success of the Quark 
Parton Model, which assumes quasi-free partons inside the hadrons. In this case 
the cross sections, if expressed in dimensionless scaling variables, are independent 
of energy. The observation of scaling in deep inelastic lepton-nucleon scattering 
led to the award of the Nobel Prize to Freedman, Kendall and Taylor in 1990. 
Even the observation of logarithmic scaling violations, both in DIS and e+e~ 
annihilation, as predicted by QCD, were observed and could be used for precise 
determinations of the strong coupling constant of QCD[^, |2l|| . 

The discovery of the beauty quark at Fermilab in Batavia(USA) in 1976 and 
the r-lepton at SLAG, both in 1976, led to the discovery of the third generation of 
quarks and leptons, of which the expected top quark is still missing. Recent LEP 
data indicate that its mass is around 166 GeV[^], thus explaining why it has not 



yet been discovered at the present accelerators. The third generation had been 
introduced into the Standard Model long before by Kobayashi and Maskawa in 
order to be able to explain the observed GP violation in the kaon system within 
the Standard Model. 
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From the total decay width of the bosons, as measured at LEP, one con- 
cludes that it couples to three different neutrinos with a mass below Mzjl ~ 45 
GeV. This strongly suggests that the number of generations of elementary par- 
ticles is not infinite, but indeed three, since the neutrinos are massless in the 



Standard Model. The three generations have been summarized in table ^TT] to- 
gether with the gauge fields, which are responsible for the low energy interactions. 

The gluons are believed to be massless, since there is no reason to assume 
that the SU (3) symmetry is broken, so one does not need Higgs fields associated 
with the low energy strong interactions. The apparent short range behaviour of 
the strong interactions is not due to the mass of the gauge bosons, but to the 
gluon self-interaction leading to confinement, as will be discussed in more detail 
afterwards. 

This chapter has been organized as follows: after a short description of the SM, 
we discuss it shortcomings and unanswered questions. They form the motivation 
for extending the SM towards a Grand Unified Theory, in which the electroweak- 
and strong forces are unified into a new force with only a single coupling constant. 
The Grand Unified Theories will be discussed in the next chapter. Although such 
unification can only happen at extremely high energies - far above the range of 
present accelerators - it still has strong implications on low energy physics, which 
can be tested at present accelerators. 



2.2 The Standard Model 

Constructing a gauge theory requires the following steps to be taken: 

• Choice of a symmetry group on the basis of the symmetry of the observed 
interactions. 

• Requirement of local gauge invariance under transformations of the sym- 
metry group. 

• Choice of the Higgs sector to introduce spontaneous symmetry breaking, 
which allows the generation of masses without breaking explicitly gauge 
invariance. Massive gauge bosons are needed to obtain the short-range be- 
haviour of the weak interactions. Adding ad-hoc mass terms, which are not 
gauge-invariant, leads to non-renormalizable field theories. In this case the 
infinities of the theory cannot be absorbed in the parameters and fields of 
the theory. With the Higgs mechanism the theory is indeed renormalizable, 
as was shown by G. 't Hooftp3|. 



Renormalization of the couplings and masses in the theory in order to relate 
the bare charges of the theory to known data. The Renormalization Group 
analysis leads to the concept of "running", i.e. energy dependent coupling 
constants, which allows the absorption of infinities in the theory into the 
coupling constants. 
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Interactions 




strong 


electro- weak 


gravitational 


unined : 


i neory 


QCD 


GSW 


quantum gravity ? 


SUGRA ? 


Symmetry 


SU{3) 


SU{2) X f/(l) 


? 


SU(5)? 


Gauge 


gi---9& 


photon 


G 


X,Y ? 


bosons 


gluons 


bosons 


graviton 


GUT bosons? 


charge 


colour 


weak isospin 
weak hypercharge 


mass 


? 



Table 2.1: The fundamental forces. The question marks indicate areas of intensive 
research. 

2.2.1 Choice of the Group Structure. 

Groups of particles observed in nature show very similar properties, thus suggest- 
ing the existence of symmetries. For example, the quarks come in three colours, 
while the weak interactions suggest the grouping of fermions into doublets. This 
leads naturally to the SU{?>) and SU{2) group structure for the strong and 
weak interactions, respectively. The electromagnetic interactions don't change 
the quantum numbers of the interacting particles, so the simple U{1) group is 
sufficient. 

Consequently, the Standard Model of the strong and electroweak interactions 
is based on the symmetry of the following unitaryQ groups: 

SU{'i)c®SU{2)L®U{l)Y. (2.1) 

The need for three colours arose in connection with the existence of hadrons 
consisting of three quarks with identical quantum numbers. According to the 
Pauli principle fermions are not allowed to be in the same state, so labeling them 
with different colours solved the problem]^. More direct experimental evidence 
for colour came from the decay width of the 7r° and the total hadronic cross 
section in e^e~ annihilation[[^. Both are proportional to the number of quark 
species and both require the number of colours to be three. 

Although colour was introduced first as an ad-hoc quantum number for the 
reasons given above, it became later evident, that its role was much more funda- 
mental, namely that it acted as the source of the field for the strong interactions 
(the "colour" field), just like the electric charge is the source of the electric field. 

The "charge" of the weak interactions is the third component of the "weak" 
isospin T3. The charged weak interactions only operate on left-handed parti- 
cles, i.e. particles with the spin aligned opposite to their momentum (negative 
helicity), so only left-handed particles are given weak isospin ±1/2 and right- 
handed particles are put into singlets (see table p.2|) . Right-handed neutrinos 
do not exist in nature, so within each generation one has 15 matter fields: 2(1) 
left(right)-handed leptons and 2x3 (2x3) left(right)-handed quarks (factor 3 for 
colour) . 

^ Unitary transformations rotate vectors, but leave their length constant. SU(N) symmetry 
groups are Special Unitary groups with determinant +1. 
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The electromagnetic interactions originate both from the exchange of the 
neutral gauge boson of the SU{2) group as well as the one from the U{1) group. 
Consequently the "charge" of the U{1) group cannot be identical with the electric 
charge, but it is the so-called weak hypercharge which is related to the 

electric charge via the Gell-Mann-Nishijima relation: 

Q = T3 + ^IV- (2.2) 

The quantum number Y\y is (B—L) for left handed doublets and 2Q for righthanded 
singlets, where the baryon number B=l/3 for quarks and for leptons, while the 
lepton number L =1 for leptons and for quarks. Since T3 and Q are conserved, 
Yw is also a conserved quantum number. The electro-weak quantum numbers 
for the elementary particle spectrum are summarized in table |2.2| . 





Generations 
1. 2. 3. 


Quantum Numbers 

Q ^3 


helicity 


L 


('■). (;). ('•). 
(;). iA (A 


1/2 -1 
-1 -1/2 -1 

2/3 1/2 1/3 
-1/3 -1/2 1/3 


R 


Gr liR Tr 

Ur Cr tR 
dR Sr bR 


-1 -2 

2/3 4/3 
-1/3 -2/3 



Table 2.2: The electro-weak quantum numbers (electric charge Q, third compo- 
nent of weak isospin T3 and weak hypercharge Yw) of the particle spectrum. The 
neutrinos u^, and u^- are the weak isospin partners of the electron(e), muon(/i) 
and tau(r) leptons, respectively. The up{u), down(c/), strange(s), charm(c), 
bottom(6) and top(t) quarks come in three colours, which have not been indi- 
cated. The primes for the left handed quarks d', s' and b' indicate the interaction 
eigenstates of the electro- weak theory, which are mixtures of the mass eigenstates, 
i.e. the real particles. The mixing matrix is the Cabibbo-Kobayshi-Maskawa ma- 
trix. The weak hypercharge Yiy equals B — L for the left-handed doublets and 
2Q for the right-handed singlets. 
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Figure 2.1: Global rotations leave the baryon colourless (a). Local rotations 
change the colour locally, thus changing the colour of the baryon (b), unless the 
colour is restored by the exchange of a gluon (c). 



















[(D®J 





Figure 2.2: Demonstration of the non-abelian character of the SU(3) rotations 
inside a colourless baryon: on the left-hand side one first exchanges a red-green 
gluon, which exchanges the colours of the quarks, and then a green-blue gluon; 
on the right-hand side the order is reversed. The final result is not the same, so 
these operations do not commute. 
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2.3 Requirement of local gauge invariance. 



The Lagrangian £ of a free fermion can be written as: 

C = ^^7/^9^^ - m^^, (2.3) 

where the first term represents the kinetic energy of the matter field \1/ with 
mass m and the second term is the energy corresponding to the mass m. The 
Euler-Lagrange equations for this C yield the Dirac equation for a free fermion. 

The unitary groups SU{N) introduced above represent rotations in N di- 
mensionalQ space. The bases for the space are provided by the eigenstates of the 
matter fields, which are the colour triplets in case of SU (3), weak isospin doublets 
in case of SU{2) and singlets for U{1). 

Arbitrary rotations of the states can be represented by 

Ar2-1 

U = exp{—id ■ F) = exp{—i ^ ■ F^) (2.4) 

k=l 

where are the rotation parameters and F^ the rotation matrices. F^ are the 
eight 3x3 Gell-Mann matrices for SU{3), denoted by A hereafter, and the well 
known Pauli matrices for SU{2) denoted by r. 

The Lagrangian is invariant under the SU{N) rotation, if £(\1'') =£(\E'), where 
\E'' = [/\E'. The mass term is clearly invariant: m\E' \E'' = m\E'f/^[/\E' = m\E'\E', 
since WU = 1 for unitary matrices. The kinetic term is only invariant under 
global transformations, i.e. transformations where is everywhere the same in 
space-time. In this case U is independent of x and can be treated as a constant 
multiplying \Ef, which leads to: '^W-f^'d^U'^ = '^WUYdf,'^ = '^^''df,'^. 

However, one could also require local instead of global gauge invariance, im- 
plying that the interactions should be invariant under rotations of the symmetry 
group for each particle separately. The motivation is simply that the interactions 
should be the same for particles belonging to the same multiplet of a symmetry 
group. For example, the interaction between a green and a blue quark should 
be the same as the interaction between a green and a red quark; therefore it 
should be allowed to perform a local colour transformation of a single quark. 
The consequence of requiring local gauge invariance is dramatic: it requires the 
introduction of intermediate gauge bosons whose quantum numbers completely 
determine the possible interactions between the matter fields, as was first shown 
by Yang and Mills in 1957 for the isopin symmetry of the strong interactions. 

Intuitively this is quite clear. Consider a hadron consisting of a colour triplet 
of quarks in a colourless groundstate. A global rotation of all quark fields will 
leave the groundstate invariant, as shown schematically in fig. p.l| . However, if a 
quark field is rotated locally, the groundstate is not colourless anymore, unless a 
"message" is mediated to the other quarks to change their colours as well. The 
"mediators" in SU{3) are the gluons, which carry a colour charge themselves and 
the local colour variation of the quark field is restored by the gluons. 

■^The SU (N) groups can be represented by x iV complex matrices A or 2N'^ real numbers. 
The unitarity requirement {A^ = A~^) imposes TV^ conditions, while requiring the determinant 
to be one imposes one more constraint, so in total the matrix is represented by iV^ — 1 real 
numbers. 
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The colour charge of the gluons is a consequence of the non-abehan character 
of SU(3), which imphes that rotations in colour space do not commute, i.e. XaXb 7^ 
AfoAa, as demonstrated in fig. p^ . If the gluons would all be colourless, they would 
not change the colour of the quarks and their exchange would be commuting. 

Mathematically, local gauge invariance is introduced by replacing the deriva- 
tive with the covariantQ derivative D^, which is required to have the following 
property: 

D'^' = [/D^, (2.5) 

i.e. the covariant derivative of the field has the same transformation properties 
as the field in contrast to the normal derivative. Clearly with this requirement 
£ is manifestly gauge invariant, since in each term of eq. p.3| the transformation 
leads to the product WU = 1 after substituting 9^ D^. 

For infinitesimal transformations the covariant derivative can be written as|^: 

D, = d,+ '-^B.Yw + ■ r + ■ A, (2.6) 

where i?^,W^ and are the field quanta ("mediators") of the U{1), SU{2) and 
S'f/(3) groups and g' , g and gs the corresponding coupling constants. 
The term • r can be explicitly written as: 

Wln + Wlr, + Wlr, = Wl(\ l]+K 




(2.7) 



The operators in the off-diagonal elements act as lowering- and raising opera- 
tors for the weak isospin. For example, they transform an electron into a neutrino 
and vice-versa, while the operator represents the neutral current interactions 
between a fermion and antifermion. 

After substituting the Gell-Mann matrices A the term G^ ■ A = Sf—iCfcAfc 
can be written similarly as: 



( Gl + - iG^ - iG^ \ 

Gl + iG^ -G' + j^G' G^'-iG'^ 



( Gl + v^a, V2G.r, \ 

V2Ggr -G^ + j=^G^ V2G,, 



\ G'^ + iG^ G' + tG' -^GS ] \ v^G,, V2G,, -^G« 

(2.8) 

This term induces transitions between the colours. For example, the off- 
diagonal element G.,.g acts like a raising operator between a green $g = (0, 1,0) 
and red $r = (1? 0, 0) field. The terms on the diagonal don't change the colour. 
Since the trace of the matrix has to be zero, there are only two independent 

^The term originates from Weyl, who tried to introduce local gauge invariance for gravity, 
thus introducing the derivative in curved space-time, which varies with the curvature, thus 
being covariant. 
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gluons, which don't change the colour. They are hnear combinations of the 
diagonal matrices A3 and Ag. 

The Wf^ gauge fields cannot represent the mediators of the weak interactions, 
since the latter have to be massive. Mass terms for W^, such as M^W^W^, 
are not gauge invariant, as can be checked from the transformation laws for the 
fields. The real fields 7, Z^, and can be obtained from the gauge fields after 
spontaneous symmetry breaking via the Higgs mechanism, as will be discussed 
in the next section. 

2.4 The Higgs mechanism. 
2.4.1 Introduction. 

The problem of mass for the fermions and weak gauge bosons can be solved by 
assuming that masses are generated dynamically through the interaction with a 
scalar field, which is assumed to be present everywhere in the vacuum, i.e. the 
space-time in which interactions take place. 

The vacuum or equivalently the groundstate, i.e. the state with the lowest 
potential energy, may have a non-zero (scalar) field value represented by $ = 
V exp{i(j)); v is called the vacuum expectation value (vev). The same minimum 
is reached for an arbitrary value of the phase 0, so there exists an infinity of 
different, but equivalent groundstates. This degeneracy of the ground state takes 
on a special significance in a quantum field theory, because the vacuum is required 
to be unique, so the phase cannot be arbitrarily at each point in space-time. Once 
a particular value of the phase is chosen, it has to remain the same everywhere, 
i.e. it cannot change locally. A scalar field with a nonzero vev therefore breaks 
local gauge invarianceQ. More details can be found in the nice introduction by 
Moriyasu0]. 

Nature has many examples of broken symmetries. Superconductivity is a well 
known example. Below the critical temperature the electrons bind into Cooper 
pairsQ The density of Cooper pairs corresponds to the vev. Owing to the weak 
binding, the effective size of a Cooper pair is large, about 10~^ cm, so every 
Cooper pair overlaps with about 10^ other Cooper pairs and this overlap "locks" 
the phases of the wave function over macroscopic distances: "Superconductivity 
is a remarkable manifestation of Quantum Mechanics on a truly macroscopic 
scale" 0. 

In the superconducting phase the photon gets an effective mass through the 
interaction with the Cooper pairs in the "vacuum", which is apparent in the 

^ An amusing analogy was proposed by A. Salam: A number of guests sitting around a 
round dinner table all have a serviette on the same side of the plate with complete symmetry. 
As soon as one guest picks up a serviette, say on the lefthand side, the symmetry is broken and 
all guests have to follow suit and take the serviette on the same side., i.e. the phases are locked 
together everywhere in the "vacuum" due to the "spontaneously broken symmetry" . 

^ The interaction of the conduction electrons with the lattice produces an attractive force. 
When the electron energies are sufficiently small, i.e. below the critical temperature, this 
attractive force overcomes the Coulomb repulsion and binds the electrons into Cooper pairs, 
in which the momenta and spins of the electrons are in opposite directions, so the Cooper pair 
forms a scalar field and its quanta have a charge two times the electron charge. 
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Figure 2.3: Shape of the Higgs potential for /i^ > (a) and /x^ < (b); (pi and 
02 are the real and imaginary parts of the Higgs field. 



Meissner effect: the magnetic field has a very short penetration depth into the 
superconductor or equivalently the photon is very massive. Before the phase 
transition the vacuum would have zero Cooper pairs, i.e. a zero vev, and the 
magnetic field can penetrate the superconductor without attenuation as expected 
for massless photons. 

This example of Quantum Mechanics and spontaneous symmetry breaking 
in superconductivity has been transferred almost literally to elementary particle 



physics by Higgs and others||2^. For the self-interaction of the Higgs field one 
considers a potential analogous to the one proposed by Ginzburg and Landau for 
superconductivity: 

V{^) = + (2.9) 

where /x^ and A are constants. The potential has a parabolic shape, if //^ > 0, but 
takes the shape of a Mexican hat for < 0, as pictured in fig. |2.3| . In the latter 
case the field free vacuum, i.e. $ = 0, corresponds to a local maximum, thus 
forming an unstable equilibrium. The groundstate corresponds to a minimum 
with a nonzero value for the field: 



1*1 = /^. (2.10) 

In superconductivity /i^ acts like the critical temperature Tc. above the 
electrons are free particles, so their phases can be rotated arbitrarily at all points 
in space, but below Tc the individual rotational freedom is lost, because the 
electrons form a coherent system, in which all phases are locked to a certain 
value. This corresponds to a single point in the minimum of the Mexican hat, 
which represents a vacuum with a nonzero vev and a well defined phase, thus 
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defining a unique vacuum. The coherent system can still be rotated as a whole 
so it is invariant under global but not under local rotations. 



2.4.2 Gauge Boson Masses and the Top Quark Mass. 

After this general introduction about the Higgs mechanism, one has to consider 
the number of Higgs fields needed to break the SU{2)l ® U{1)y symmetry to 
the U{l)em symmetry. The latter must have one massless gauge boson, while the 
W and Z bosons must be massive. This can be achieved by choosing $ to be a 
complex SU{2) doublet with definite hypercharge {Yw = 

In order to understand the interactions of the Higgs field with other particles, 
one considers the following Lagrangian for a scalar field: 

Ch = {D^^)^{D^^)-V{^). (2.12) 

The first term is the usual kinetic energy term for a scalar particle, for which the 
Euler- Lagrange equations lead to the Klein-Gordon equation of motion. Instead 



of the normal derivative, the covariant derivative is used in eq. |2.12| in order to 
ensure local gauge invariance under SU{2) ® U{1) rotations. 

The vacuum is known to be neutral. Therefore the groundstate of $ has to 
be of the form (0,v). Furthermore $(x) has to be constant everywhere in order 
to have zero kinetic energy, i.e. the derivative term in Ch disappears. 

The quantum fluctuations of the field around the ground state can be parametrised 
as follows, if we include an arbitrary SU (2) phase factor: 



* = e'^^^"^" ut N • (2.13) 




The (real) fields C,{x) are excitations of the field along the potential minimum. 
They correspond to the massless Goldstone bosons of a global symmetry, in this 
case three for the three rotations of the SU (2) group. However, in a local gauge 
theory these massless bosons can be eliminated by a suitable rotation: 



= e-<^-yf^^^) = { /; . (2.14) 




Consequently the field C, has no physical significance. Only the real field h{x) can 
be interpreted as a real (Higgs) particle. The original field $ with four degrees of 
freedom has lost three degrees of freedom; these are recovered as the longitudinal 
polarizations of the three heavy gauge bosons. 

The kinetic part of eq. |2.12| gives rise to mass terms for the vector bosons, 
which can be written as {Yy/ = 1): 



C 



(2.15) 



H 
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or substituting for $ its vacuum expectation value v one obtains from the off- 
diagonal terms (by writing the r matrices explicitly, see eq. p.7|) 

(f)'KT + (W^.')') (2-16) 
and from the diagonal terms: 

(2.17) 

Since mass terms of physical fields have to be diagonal, one obtains the "physical" 
gauge fields of the broken symmetry by diagonalizing the mass term: 

{BlW'^^) U-'U M U-'U ^ (2.18) 

where U represents a unitary matrix 

^ 1 / 9 9' \ ^ ( cos 9w sinOw \ -|q\ 
VF^T7 [-9' 9 )~[- sin^H/ cosOw J ^ ' 

Consequently the real fields become a mixture of the gauge fields: 

( ) = ^ ( ) P-^") 

and the matrix UMU~^ becomes a diagonal matrix for a suitable mixing angle 
In these fields the mass terms have the form 

M'^W;W-^ + \{A,^Z,)^l Ml)[%) (2-21) 

with 

= \g^v^ (2.22) 

M| = ^ v\ (2.23) 

Here v is the vacuum expectation value of the Higgs potential, which for the 
known gauge boson masses and couplings can be calculated to be[|: 

V ^ 174 GeV. (2.24) 

The neutral part of the Lagrangian, if expressed in terms of the physical fields, 
can be written as: 

CZt'" = -M9' cos OwieR-f^eR + ^Vlj^ul + ^CLl^eL) - sinOwiViYuL - eLl^CL)] 

+Z^,[g' sm9w{eRj^eR + ]^Ll^yL + ^^lI^gl) + ]^g cos 9 wi^Ll'^i^L - eLT^e^)] 

(2.25) 



^Sometimes the Higgs field is normalized by 1/V2, in which case v ~ 256 GeV. 
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Figure 2.4: Geometric picture of the relations between the electroweak couphng 
constants. 



The photon field should only couple to the electron fields and not to the neutrinos, 
so the terms proportional to g' cos 6w and g sin 6w should cancel and the coupling 
to the electrons has to be the electric charge e. This can be achieved by requiring: 

g' cosOw = g sin 9w = G . (2.26) 

Hence 

g^'^ gg^ 

tan9w = — ; sin^ 9w = — ^ and e = —====. (2.27) 
9 9^ + 9^ V9 + 9 

A geometric picture of these relations is shown in fig. |2.4| . From these rela- 
tions and the relations between masses and couplings ( |2.22| and |2.23| ) one finds 
the famous relation between the electroweak mixing angle and the gauge boson 
masses: 

Mw = cosOw ■ Mz or sin^ = I - ^ ■ (2.28) 

The value of M\y can also be related to the precisely measured muon decay 
constant = 1.16639(2) ■ 10"^ GeV-^. If calculated in the SM, one finds: 

71 " 8sin2 OwM^' 

This relation can be used to calculate the gauge boson masses from measured 
coupling constants a, G^ and smOw'- 

o vra 1 
Mw = 27r- 2.30 

^ V^Gf, sin^ 9w cos^ 9w ^ ^ 

(2.32) 

Inserting sin^ 9w = 0.23 and 1/a = 137.036 yields Mz=SS GeV. However, these 
relations are only at tree level. Radiative corrections depend on the as yet un- 
known top mass. Fitting the unknown top mass to the measured Mz mass, the 
electroweak asymmetries and the cross sections at LEP yields p2[| : 

Mtop = 166;!;^g (stat.) 1^11 {unknown Higgs). (2.33) 
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Also the fermions can interact with the scalar field, albeit not necessarily with 
the gauge coupling constant. The Lagrangian for the interaction of the leptons 
with the Higgs field can be written as: 



Substituting the vacuum expectation value for $ yields 

eR + cr (0, v) 



V2 









C) 








r_ , - 1 

[eieR + eRCL] 



V2 



(2.34) 



-ee 



V2 
(2.35) 

The Yukawa coupling constant gy is a free parameter, which has to be adjusted 
such that rrif. = gyVy/^. Thus the coupling gy is proportional to the mass of 
the particle and consequently the coupling of the Higgs field to fermions is pro- 
portional to the mass of the fermion, a prediction of utmost importance to the 
experimental search for the Higgs boson. 

Note that the neutrino stays massless with the choice of the Lagrangian, since 
no mass term for the neutrino appears in eq. p.35| . 



2.4.3 Summary on the Higgs mechanism. 

In summary, the Higgs mechanism assumed the existence of a scalar field $ = 
$0 exp {i6{x)). After spontaneous symmetry breaking the phases are "locked" 
over macroscopic distances, so the field averaged over all phases is not zero any- 
more and $ develops a vacuum expectation value. The interaction of the fermions 
and gauge bosons with this coherent system of scalar fields $ gives rise to effec- 
tive particle masses, just like the interaction of the electromagnetic field with the 
Cooper pairs inside a superconductor can be described by an effective photon 
mass. 

The vacuum corresponds to the groundstate with minimal potential energy 
and zero kinetic energy. At high enough temperatures the thermal fluctuations of 
the Higgs particles about the groundstate become so strong that the coherence is 
lost, i.e. $(x) = constant is not true anymore. In other words a phase transition 
from the ground state with broken symmetry ($ ^ 0) to the symmetric ground- 
state takes place. In the symmetric phase the groundstate is invariant again under 
local SU{2) rotations, since the phases can be adjusted locally without changing 
the groundstate with < $ >= 0. In the latter case all masses disappear, since 
they are proportional to < $ >= 0. 

Both, the fermion and gauge boson masses are generated through the interac- 
tion with the Higgs fleld. Since the interactions are proportional to the coupling 
constants, one flnds a relation between masses and coupling constants. For the 
fermions the Yukawa coupling constant is proportional to the fermion mass and 
the mass ratio of the W and Z bosons is only dependent on the electroweak 
mixing angle (see eq. p.28|) . This mass relation is in excellent agreement with ex- 



perimental data after including radiative corrections. Hence, it is the first indirect 
evidence that the gauge bosons masses are indeed generated by the interaction 
with a scalar field, since otherwise there is no reason to expect the masses of the 
charged and neutral gauge bosons to be related in such a specific way via the 
couplings. 



18 



QED QCD 

Figure 2.5: The effective charge distribution around an electric charge (QED) and 
colour charge (QCD). At higher one probes smaller distances, thus observing 
a larger (smaller) effective charge, i.e. a larger (smaller) coupling constant in 
QED (QCD). 
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Figure 2.6: Running of the three couphng constants in the Standard Model owing 
to the different space charge distributions (compare fig. 2.5. 

2.5 Running Coupling Constants 

In a Quantum Field Theory the couphng constants are only effective constants 
at a certain energy. They arc energy, or equivalently distance dependent through 
virtual corrections, both in QED and in QCD. 

However, in QED the coupling constant increases as function of Q^, while in 
QCD the coupling constant decreases. A simple picture for this behaviour is the 
following: 

• The electric field around a pointlike electric charge diverges like 1/r. In 
such a strong field electron-positron pairs can be created with a lifetime 
determined by Heisenberg's uncertainty relations. These virtual e^e~ pairs 
orient themselves in the electric field, thus giving rise to vacuum polariza- 
tion, just like the atoms in a dielectric are polarized by an external electric 
field. This vacuum polarization screens the "bare" charge, so at a large 
distance one observes only an effective charge. This causes deviations from 
Coulomb's law, as observed in the well-known Lamb shift of the energy 
levels of the hydrogen atom. If the electric charge is probed at higher en- 
ergies (or shorter distances), one penetrates the shielding from the vacuum 
polarization deeper and observes more of the bare charge, or equivalently 
one observes a larger coupling constant. 

• In QCD the situation is more complicated: the colour charge is surrounded 
by a cloud of gluons and virtual qq pairs; since the gluons themselves 
carry a colour charge, one has two contributions: a shielding of the bare 
charge by the qq pairs and an increase of the colour charge by the gluon 
cloud. The net effect of the vacuum polarization is an increase of the total 
colour charge, provided not too many qq pairs contribute, which is the case 
if the number of generations is below 16, see hereafter. If one probes this 
charge at smaller distances, one penetrates part of the "antishielding" , thus 
observing a smaller colour charge at higher energies. So it is the fact that 
gluons carry colour themselves which makes the coupling decrease at small 



20 




screening antiscreening 




Figure 2.7: Loop corrections in QED (a) and QCD (b). In QED only the fermions 
contribute in the loops, which causes a screening of the bare charge. In QCD 
also the bosons contribute through the gluon selfinteraction, which enhances the 
bare charge. This antiscreening dominates over the screening. 



distances (or high energies). This property is called asymptotic freedom 
and it explains why in deep inelastic lepton-nucleus scattering experiments 
the quarks inside a nucleus appear quasi free in spite of the fact that they 
are tightly bound inside a nucleus. The increase of at large distances 
explains qualitatively why it is so difficult to separate the quarks inside a 
hadron: the larger the distance the more energy one needs to separate them 
even further. If the energy of the colour field is too high, it is transformed 
into mass, thus generating new quarks, which then recombine with the old 
ones to form new hadrons, so one always ends up with a system of hadrons 
instead of free quarks. 

The space charges from the virtual pairs surrounding an electric charge and 
colour charge are shown schematically in fig. The different vacuum polar- 

izations lead to the energy dependence of the coupling constants sketched in fig. 
|2.6| . The colour field becomes infinitely dense at the QCD scale A ^ 200 MeV 
(see hereafter). So the confinement radius of typical hadrons is (9(200 MeV) or 
one Fermi (10~^^ cm). 

The vacuum polarization effects can be calculated from the loop diagrams to 
the gauge bosons. The main difference between the charge distribution in QED 
and QCD originates from the diagrams shown in fig. |2.7|. In addition one has 



to consider diagrams of the type shown in fig. p.8| . The ultraviolet divergences 
(Q^ oo) in these diagrams can be absorbed in the coupling constants in a 
renormalizable theory. All other divergences are canceled at the amplitude level 
by summing the appropriate amplitudes. The first step in such calculations is 
the regularization of the divergences, i.e. separating the divergent parts in the 
mathematical expressions. The second step is the renormalization of physical 
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Figure 2.8: First order vacuum polarization diagrams. 

quantities, like charge and mass, to absorb the divergent parts of the amplitudes, 
i.e. replace the "bare" quantities of the theory with measured quantities. For 
example, the loop corrections to the photon propagator diverge, if the momentum 
transfer k in the loop is integrated to infinity. If one introduces a cutoff /xq for 
large values of k, one finds for the regularized amplitude of the sum of the Born 



term Mq and the loop corrections] 26 



M' = e'(l-—\n^] fl + ^ln^^ for » (2.36) 
y 377 m'^ ) \ 377 m^y 

The divergent part depending on the cutoff parameter /xq disappears, if one re- 
places the "bare" charge e by the renormalized charge cr. 

(2.37) 

i.e. the "bare" charge, occurring in the Dirac equation, is renormalized to a 
measurable quantity cr. For one usually takes the Thomson limit for Compton 
scattering, i.e. 76 — > 76 for k 

„2 




6T( a , , 

= — — (2.38) 
3 mi 



e 

^2 



with a = eyAn = 1/137.036 and = 0.00051 GeV. 

After regularization and renormalization to a measured quantity (in this case 
using the so called "on shell" scheme, i.e. one uses the mass and charge of a 
free electron as measured at low energy), one is left with a dependent but 
finite part of the vacuum polarization. This can be absorbed in a dependent 
coupling constant, which in case of QED becomes for » m?: 

a{Q^) = a(l + ^\n^\ (2.39) 
V Sit mi 
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If one sums more loops, this yields terms (^)"(ln^)™' and retaining only the 
leading logarithms (i.e. n=m), these terms can be summed to: 



a 



(1-f In^ 



(2.40) 



smce 



Ex" 

n=0 



1 — X 



(2.41) 



Of course, the total dependence is obtained by summing over all possible 
fermion loops in the photon propagator. 

These vacuum polarization effects are non- negligible. For example, at LEP 
accelerator energies a has increased from its low energy value 1/137 to 1/128 or 
about 6%. 

The diagrams of fig. p.7| b yield similarly to eq. p.40| : 



1 + 



47r 



11 



/i2 



(2.42) 



Note that decreases with increasing if 11 — 2Nf/3 > or Nf < 16, 
thus leading to asymptotic freedom at high energy. This is in contrast to the 
dependence of a(Q^) in eq. p.40| , which increases with increasing Q^. Since 
as becomes infinite at small Q^, one cannot take this scale as a reference scale. 
Instead one could choose as renormalization point the "confinement scale" A, i.e. 
as — > oo, if A. In this case eq. 2.42 becomes independent of /i, since the 1 

in brackets becomes negligible, so one obtains: 



Air 



'11 



(2.43) 



The definition of A depends on the renormalization scheme. The most widely 



used scheme is the MS scheme p7|] , which we will use here. Other schemes can 



be used as well and simple relations between the definitions of A exist] 28 



The higher order corrections are usually calculated with the renormalization 
group technique, which yields for the /i dependence of a coupling constant a : 



Qa 



(2.44) 



The first two terms in this perturbative expansion are renormalization-scheme 
independent. Their specific values are given in the appendix. The first order 
solution of eq. |2.44| is simple: 



a(g^ 



1 1 (Q'^ 



(2.45) 



where Qq is a reference energy. One observes a linear relation between the change 
in the inverse of the coupling constant and the logarithm of the energy. The slope 
depends on the sign of /3o, which is positive for QED, but negative for QCD, thus 
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leading to asymptotic freedom in the latter case. The second order corrections 
are so small, that they do not change this conclusion. Higher order terms depend 
on the renormalization prescription. In higher orders there are also corrections 
from Higgs particles and gauge bosons in the loops. Therefore the running of 
a given coupling constant depends slightly on the value of the other coupling 
constants and the Yukawa couplings. These higher order corrections cause the 
RGE equations to be coupled, so one has to solve a large number of coupled 
differential equations. All these equations are summarized in the appendix. 
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Chapter 3 

Grand Unified Theories. 



3.1 Motivation 

The Standard Model describes all observed interactions between elementary par- 
ticles with astonishing precision. Nevertheless, it cannot be considered to be 
the ultimate theory because the many unanswered questions remain a problem. 
Among them: 

• The Gauge Problem 

Why are there three independent symmetry groups? 

• The Pcirameter Problem 

How can one reduce the number of free parameters? (At least 18 from 
the couplings, the mixing parameters, the Yukawa couplings and the Higgs 
potential.) 

• The Fermion Problem 

Why are there three generations of quarks and leptons? What is the origin 
of the the symmetry between quarks and leptons? Are they composite 
particles of more fundamental objects? 

• The Charge Quantization Problem 

Why do protons and electrons have exactly opposite electric charges? 

• The Hierarchy Problem 

Why is the weak scale so small compared with the GUT scale, i.e. why is 

• The Fine-tuning Problem 

Radiative corrections to the Higgs masses and gauge boson masses have 
quadratic divergences. For example, AM^ ^ ^i^pianck)- other words, 
the corrections to the Higgs masses are many orders of magnitude larger 
than the masses themselves, since they arc expected to be of the order of 
the electroweak gauge boson masses. This requires extremely unnatural 
fine-tuning in the parameters of the Higgs potential. This "fine-tuning" 
problem is solved in the supersymmetric extension of the SM, as will be 
discussed afterwards. 
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3.2 Grand Unification 



The problems mentioned above can be partly solved by assuming the symmetry 
groups SU{3)c ® SU{2)l ® U{X)y are part of a larger group G, i.e. 



G D SU{'i)c ® SU{2)l ® U{1)y. 



(3.1) 



The smallest group G is the SU(5) groupflP^I , so the minimal extension of the 
SM towards a GUT is based on the SU{b) group. Throughout this paper we will 
only consider this minimal extension. The group G has a single coupling constant 
for all interactions and the observed differences in the couplings at low energy are 
caused by radiative corrections. As discussed before, the strong coupling constant 
decreases with increasing energy, while the electromagnetic one increases with 
energy, so that at some high energy they will become equal. Since the changes 
with energy are only logarithmic (eq. p. 451 ), the unification scale is high, namely 
of the order of 10^^ — 10^^ GeV, depending on the assumed particle content in 
the loop diagrams. 

In the SU (5) group the 15 particles and antiparticles of the first generation 
can be fit into the 5-plet0 and 10-plet: 
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(3.2) 



+ Ug 

V +dg 

The superscript G indicates the charge conjugated particle, i.e. the antiparticle 
and all particles are chosen to be left-handed, since a left-handed antiparticle 
transforms like a right-handed particle. Thus the superscript G implies a right- 
handed singlet with weak isospin equal zero. 

With this multiplet structure the sum of the quantum numbers Q, T3 and Y 
is zero within one multiplet, as required, since the corresponding operators are 
represented by traceless matrices. 

Note that there is no space for the antineutrino in these multiplets, so within 
the minimal SU{5) the neutrino must be massless, since for a massive particle 
the right-handed helicity state is also present. Of course, it is possible to put a 
right-handed neutrino into a singlet representation. 

SU{5) rotations can be represented by 5 x 5 matrices. Local gauge invariance 
requires the introduction of iV^ — 1 = 24 gauge fields (the "mediators"), which 
cause the interactions between the matter fields. The gauge fields transform 
under the adjoint representation of the SU{5) group, which can be written in 



cannot be the direct product of the SU{3), SU{2) and U{1) groups, since this would not 
represent a new unified force with a single coupling constant, but still require three independent 
coupling constants. 

^The bar indicates the complementary representation of the fundamental representation. 



26 



Figure 3.1: GUT proton decays through the exchange of X and Y gauge bosons. 



matrix form as (compare eqns. |2.7| and 
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(3.3) 



The G's represent the gluon fields of equation p.8| , while the W^s and -B's are the 
gauge fields of the SU{2) symmetry groups. The X and F's are new gauge bosons, 
which represent interactions, in which quarks are transformed into leptons and 
vice- versa, as should be apparent if one operates with this matrix on the 5-plet. 
Consequently, the X {Y) bosons, which couple to the electron (neutrino) and 
(i-quark must have electric charge 4/3 (1/3). 



3.3 5(7(5) predictions 
3.3.1 Proton decay 

The X and Y gauge bosons can introduce transitions between quarks and leptons, 
thus violating lepton and baryon number^. This can lead to the following proton 
and neutron decays (see fig. p.l|) : 

p — s> e^7r° n — e+vr^ 

p e^p^ n — s> e+p^ 

p e^uo*^ n — s> vuo^ 

p e^rj n Vti^ (3.4) 

p — s> Vtt^ n 



p up^ 
p -> Vf,K+ 



The decays with kaons in the final state are allowed through flavour mixing, 
i.e. the interaction eigenstates are not necessarily the mass eigenstates. 



^The difference between lepton and baryon number B-L is conserved in these transitions. 
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For the lifetime of the nucleon one writes in analogy to muon decay: 



The proton mass irip to the fifth power originates from the phase space in case the 
final states are much lighter than the proton, which is the case for the dominant 
decay mode: p —>■ e~^7iQ. After this prediction of an unstable proton in grand 
unified theories, a great deal of activity developed and the lower limit on the 
proton life time increased toPD| 



Tp > 5 ■ lO'^^ yrs (3.6) 



for the dominant decay mode p e+vr". From equation p.5| this implies (for 
as = 1/24, see chapter |^) 

Mx > 10^^ GeV. (3.7) 

From the extrapolation of the couplings in the SU (5) model to high energies 
one expects the unification scale to be reached well below 10^^ GeV, so the proton 
lifetime measurements exclude the minimal SU{5) model as a viable GUT. As 
will be discussed later, the supersymmetric extension of the SU{5) model has the 
unification point well above 10^^ GeV. 

3.3.2 Baryon Asymmetry 

The heavy gauge bosons responsible for the unified force cannot be produced 
with conventional accelerators, but energies above 10^^ were easily accessible 
during the birth of our universe. This could have led to an excess of matter over 
antimatter right at the beginning, since the X and Y bosons can decay into pure 
matter, e.g. X — > uu, which is allowed because the charge of the X boson is 
4/3. As pointed out by Sakharov [pI| such an excess is possible if both C and 
CP are violated, if the baryon number B is violated, and if the process goes 
through a phase of non-equilibrium. All three conditions are possible within the 
SU (5) model. The non-equilibrium phase happens if the hot universe cools down 
and arrives at a temperature, too low to generate X and Y bosons anymore, so 
only the decays are possible. Since the CP violation is expected to be small, the 
excess of matter over antimatter will be small, so most of the matter annihilated 
with antimatter into enormous number of photons. This would explain why the 
number of photons over baryons is so large: 

^ ^ 10^° (3.8) 

However, later it was realized that the electroweak phase transition may wash 
out any (B-|-L) excess generated by GUT's. One then has to explain the observed 



baryon asymmetry by the electroweak baryogenesis, which is actively studied ||32[| 
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Figure 3.2: Radiative corrections to particle masses. 



3.3.3 Charge Quantization 

From the fact that quarks and leptons are assigned to the same multiplet the 
charges must be related, since the trace of any generator has to be zero. For 
example, the charge operator Q on the fundamental representation yields: 

TrQ = Tr{q^,q^,q^,e,0) = (3.9) 

or in other words, in SU{5) the electric charge of the d-quark has to be 1/3 of 
the charge of an electron! Similarly, one finds the charge of the ti-quark is 2/3 of 
the positron charge, so the total charge of the proton (=uud) has to be exactly 
opposite to the charge of an electron. 



3.3.4 Prediction of sin^6' 



w 



If the SU{2) and U{1) groups have equal coupling constants, the electroweak 
mixing angle can be calculated easily, since it is given by the ratio g' ^/ {g'^ + g' ^) 



(see eq. |2.27|) , which would be 1/2 for equal coupling constants. However, the 
argument is slightly more subtle, since for unitary transformations the rotation 
matrices have to be normalized such that 

TrFkFi = 6ki. (3.10) 

This normalization is not critical in case one has independent coupling constants 
for the subgroups, since a "wrong" normalization for a rotation matrix can al- 
ways be corrected by a redefinition of the corresponding coupling constants, as 
is apparent from equation This freedom is lost, if one has a single coupling 



constant, so one has to be careful about the relative normalization. It turns out, 
that the Gell-Mann and Pauli rotation matrices of the SU{3) and SU{2) groups 
have the correct normalization, but the normalization of the weak hypercharge 
operator needs to be changed. Defining l/2Yw = CTq and substituting this into 
the Gell-Mann-Nishijima relation |2.2| yields: 



Q = Ts + CTo (3.11) 
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Requiring the same normalization for T3 and Tq implies from equation |3.10| : 

Tr Q"^ = {1 + C^) Tr (3.12) 
or inserting numbers from the 5-plet of SU (5) yields: 



,2 TrQ2 3-1/9 + 1 8 



Replacing in the covariant derivative (eq. p.6|) l/2Yw with CTq implies 
g'CTo = g^To or: 

95 = Cg', (3.14) 



where = 5/3 from eq. 3.13 . With this normalization the electroweak mixing 



angle after unification becomes: 

The manifest disagreement with the experimental value of 0.23 at low energies 
brought the SU{^) model originally into discredit, until it was noticed that the 
running of the couplings between the unification scale and low energies could 
reduce the value of sin^ 6w considerably. As we will show in the last chapter, 
with the very precise measurement of sin^ 6^ at LEP, unification of the three 
coupling constants within the SU{b) model is excluded, and just as in the case 
of the proton life time, supersymmetry comes to the rescue and unification is 
perfectly possible within the supersymmetric extension of SU{b). 

Note that the prediction of sin^ = 3/8 is not specific to the SU{b) model, 
but is true for any group with SU{i)c ® SU{2)l ® U{1)y as subgroups, implying 
that Q, T3 and Yy/ are generators with traces equal zero and thus leading to the 
predictions given above. 



3.4 Spontaneous Symmetry Breaking in SU{5) 



The SU{5) symmetry is certainly broken, since the new force corresponding to 
the exchange of the X and Y bosons would lead to very rapid proton decay, if 
these new gauge bosons were massless. As mentioned above, from the limit on 
the proton life time these SU (5) gauge boson have to be very heavy, i.e. masses 
above 10^^ GeV. The generation of masses can be obtained again in a gauge 
invariant way via the Higgs mechanism. The Higgs field is chosen in the adjoint 
representation 24 and the minimum < $24 > can be chosen in the following way: 



/ 1 



< $24 >= V24 



\ 



\ 



(3.16) 



The 12 X,Y gauge bosons of the SU (5) group require a mass: 

25 



Mi 



M2 



2 2 

-95^24 



(3.17) 
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after 'eating' 12 of the 24 scalar fields in the adjoint representation, thus providing 
the longitudinal degrees of freedom. The field $24 is invariant under the rotations 
of the SU{3)c ® SU{2)l (S> U{1)y group, so this symmetry is not broken and the 
corresponding gauge bosons, including the W and Z bosons, remain massless. 
after the first stage of SU (5) symmetry breaking. 

The usual breakdown of the electroweak symmetry to SU{3)c ® f^(l)em is 
achieved by a 5-plet $5 of Higgs fields, for which the minimum of the effective 
potential can be chosen at: 



< $5 >= V5 



( \ 




VI/ 



(3.18) 



The fourth and fifth component of $5 correspond to the SU{2) doublet 
of the SM (see eq. |2.11| ). Since the total charge in a representation has to be zero 
again, the first triplet of complex fields in $5, which transforms as (3,l)_2/3 and 
(3*, 1)2/3, must have charge |1/3|. Since they couple to all fermions with mass, 
they can induce proton decay: 



+ u 
Vf. + d 



(3.19) 



Such decays can be suppressed only by sufficiently high masses of the coloured 
Higgs triplet. These can obtain high masses through interaction terms between 
$5 and $24. 

Note that from eq. |3.17| < $24 > has to be of the order of Mx, while $5 has 
to be of the order of Mw, since 



and 



w 



w 



cos 9 



w 



(3.20) 



(3.21) 



or more precisely = 1/ ^Gf= 174 GeV. The minimum of the Higgs potential 
involves both $5 and $24- Despite this mixing, the ratio f5/f24 ~ 10^^^ has to 
be preserved (hierarchy problem). Radiative corrections spoil usually such a fine- 
tuning, so SU{J)) is in trouble. As will be discussed later, also here supersymmetry 
offers solutions for both this fine-tuning and the hierarchy problem. 



3.5 Relations between Quark and Lepton Masses 

The Higgs 5-plet $5 can be used to generate fermion masses. Since the 5-plet of 
the matter fields contains both leptons and down-type quarks, their masses are 
related, while the up-type quark masses are free parameters. At the GUT scale 
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one expects: 



rrid 
rus 
rUb 



(3.22) 
(3.23) 
(3.24) 



Unfortunately the masses of the hght quarks have large uncertainties from the 
binding energies in the hadrons, but the b-quark mass can be correctly predicted 
from the r-mass after including radiative corrections (see fig. ^.2| for typical 
graphs) . 

Since the corrections from graphs involving the strong coupling constant 



are dominant, one expects in first order []33| 



O 



a 



0{3) 



(3.25) 



More precise formulae are given in the appendix and will be used in the last 
chapter in a quantitative analysis, since the b-quark mass gives a rather strong 
constraint on the evolution of the couplings and through the radiative corrections 
involving the Yukawa couplings on the top quark mass. 
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Chapter 4 
Supersymmetry 



4.1 Motivation 

Supersymmetry 1 13, 14 1 presupposes a symmetry between fermions and bosons, 
which can be reahzed in nature only if one assumes each particle with spin j 
has a supersymmetric partner with spin j-1/2. This leads to a doubling of the 



particle spectrum (see table |4.1|), which are assigned to two supermultiplets: the 
vector multiplet for the gauge bosons and the chiral multiplet for the matter 
fields. Unfortunately the supersymmetric particles or "sparticles" have not been 
observed so far, so either supersymmetry is an elegant idea, which has nothing 
to do with reality, or supersymmetry is not an exact symmetry, in which case the 
sparticles can be heavier than the particles. Many people opt for the latter way 
out, since there are many good reasons to believe in supersymmetry: 

• SUSY solves the fine-tuning problem 



As mentioned before, the radiative corrections in the SU{5) model have 
quadratic divergences from the diagrams in fig. |]8|, which lead to AM^ ^ 
0{Mx), where Mx is a cutoff scale, typically the unification scale if no 
other scales introduce new physics beforehand. 

However, in SUSY the loop corrections contain both fermions (F) and 
bosons (B) in the loops, which according to the Feynman rules contribute 



VECTOR MULTIPLET 


CHIRAL MULTIPLET 


J= 1 


J =1/2 


J= 1/2 


J = 


9 
B 


9 
B 


Hi,H2 


Hi, H2 



Table 4.1: Assignment of gauge fields to the vector superfield and the matter 
fields to the chiral superfield. 
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with an opposite sign, i.e. 



AMfj ^ 0{a) \Ml - Ml\ ^ O{10~^)MIusy (4-1) 

where Msusy is a typical SUSY mass scale. In other words, the fine-tuning 
problem disappears, if the SUSY partners are not too heavy compared with 
the known fermions. An estimate of the required SUSY breaking scale can 
be obtained by considering that the masses of the weak gauge bosons and 
Higgs masses are both obtained by multiplying the vacuum expectation 
value of the Higgs field (see previous chapter) with a coupling constant, so 
one expects Mw ~ Mh- Requiring that the radiative corrections are not 
much larger than the masses themselves, i.e. AMy/ < Myy, or replacing 
Mw by Mh, AMh < C(102), yields after substitution into eq. ^ 

Msusy < 10^ GeV. (4.2) 



SUSY offers a solution for the hierarchy problem 

The possible explanation for the small ratio M^/M"^ ^ 10~^^ is simple 
in SUSY models: large radiative corrections from the top-quark Yukawa 
coupling to the Higgs sector drive one of the Higgs masses squared nega- 
tive, thus changing the shape of the effective potential from the parabolic 
shape to the Mexican hat (see fig. ^.4| ) and triggering electroweak symme- 
try breaking|3^. Since radiative corrections are logarithmic in energy, this 
automatically leads to a large hierarchy between the scales. In the SM one 
could invoke a similar mechanism for the triggering of electroweak symme- 
try breaking, but in that case the quadratic divergences in the radiative 
corrections would upset the argument. 

In the MSSM the electroweak scale is governed by the starting values of 
the parameters at the GUT scale and the top-quark mass. This strongly 
constrains the SUSY mass spectrum, as will be discussed in the last chapter. 

SUSY yields unification of the coupling constants 

After the precise measurements of the SU{3)c ® SU{2)l ® U{1)y coupling 
constants, the possibility of coupling constant unification within the SM 
could be excluded, since after extrapolation to high energies the three cou- 
pling constants would not meet in a single point. This is demonstrated in 



the upper part of fig. which shows the evolution of the inverse of the 
couplings as function of the logarithm of the energy. In this presentation 
the evolution becomes a straight line in first order, as is apparent from the 
solution of the ROE (eqn. |2.45| ). The second order corrections, which have 
been included in fig. |4.1| by using eqs. [A.20| from the appendix, are so 



small, that they cause no visible deviation from a straight line. 

A single unification point is excluded by more than 8 standard deviations. 
The curve l/as meets the crossing point of the other two coupling constants 
only for a starting value at as{Mz) = 0.07, while the measured value is 



0.12 ± 0.006|22|. This is an exciting result, since it means unification can 
only be obtained, if new physics enters between the electroweak and the 
Planck scale! 



34 



Unification of the Couplings of the 
Electromagnetic, Weak and Strong Forces 




0'^ liiGey] 



t 1 10 
1 0' 1 0'GeV 



Minimal 
Supersymmetric 
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0' 



10 t ^[GeV] 



Figure 4.1: Evolution of the inverse of the three couphng constants in the Stan- 
dard Model (SM) (top) and in the supersymmetric extension of the SM (MSSM) 
(bottom). Only in the latter case unification is obtained. The SUSY particles are 
assumed to contribute only above the effective SUSY scale Msusy of about one 
TeV, which causes the change in slope in the evolution of the couplings. The 68% 
C.L. for this scale is indicated by the vertical hues (dashed). The evolution of the 
couplings was calculated in second order (see section A. 2 of the appendix with 
the constants (3i and j3ij calculated for the MSSM above Msusy in the bottom 
part and for the SM elsewhere). The thickness of the lines represents the error 
in the coupling constants. oc 
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MsusY [GeV] Mgut IGg\/1 



Figure 4.2: distribution for Msusy a-iadMcur- 



It turns out that within the SUSY model perfect unification can be obtained 
if the SUSY masses are of the order of one TeV. This is shown in the bottom 



part of fig. PI the SUSY particles are assumed to contribute effectively 
to the running of the coupling constants only for energies above the typical 
SUSY mass scale, which causes the change in the slope of the lines near one 
TeV. From a fit requiring unification one finds for the breakpoint Msusy 



and the unification point Mc^rpSl, p6 



Msusy = 10^ ' ^ ^ GeV (4.3) 
Mgut = lO^^-^ ^ ^ GeV (4.4) 
aauT'^ = 26.3 ± 1.9 ± 1.0, (4.5) 

where acur = oi/'^'^- The first error originates from the uncertainty in 
the coupling constant, while the second error is due to the uncertainty in 
the mass splittings between the SUSY particles. The distributions of 



Msusy and Mgut for the fit in the bottom part of fig. [4.1| are shown in fig. 

These figures are an update of the published figures using the newest 
values of the coupling constants, as shown in the figure pBj . 

Note that the parametrisation of the SUSY mass spectrum with a single 
mass scale is not adequate and leads to uncertainties. However, the errors 
in the coupling constants (mainly in a^) are large and the uncertainties from 
mass splittings between the sparticles are more than a factor two smaller 
(see eq. |4.3| ). In the last chapter the unification including a more detailed 
treatment of the mass splittings will be studied. 

One can ask: 'What is the significance of this observation? For many people 
it was the first "evidence" for supersymmetry, especially since Msusy was 
found in the range where the fine-tuning problem does not reappear (see eq. 
Consequently the results triggered a revival of the interest in SUSY, 
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as was apparent from the fact that ref. |TT|] with the fit of the unification 



of the couphng constants, as exemphfied in figs. |4.1| and reached the 



Top- Ten of the citation hst, thus leading to discussions in practically all 



popular journals|37 



Non-SUSY enthusiasts were considering unification obvious: with a total 
of three free parameters {Mqut, c^gut and Msusy) and three equations 
one can naively always find a solution. The latter statement is certainly 
not true: searching for other types of new physics with the masses as free 
particles yields only rarely unification, especially if one requires in addition 
that the unification scale is above 10^^ GeV in order to be consistent with 
the proton lifetime limits and below the Planck scale in order to be in the 
regime where gravity can be neglected. From the 1600 models tried, only 
a handful yielded unification!^. The reason is simple: introducing new 



particles usually alters all three couplings simultaneously, thus giving rise 
to strong correlations between the slopes of the three lines. For example, 
adding a fourth family of particles with an arbitrary mass will never yield 
unification, since it changes the slopes of all three coupling by the same 
amount, so if with three families unification cannot be obtained, it will not 
work with four families either, even if one has an additional free parameter! 
Nevertheless, unification does not prove supersymmetry, it only gives an 
interesting hint. The real proof would be the observation of the sparticles. 

Unification with gravity 



The space-time symmetry group is the Poincare group. Requiring local 
gauge invariance under the transformations of this group leads to the Ein- 
stein theory of gravitation. Localizing both the internal and the space-time 
symmetry groups yields the Yang-Mills gauge fields and the gravitational 
fields. This paves the way for the unification of gravity with the strong 
and electroweak interactions. The only non-trivial unification of an inter- 
nal symmetry and the space-time symmetry group is the supersymmetry 
group, so supersymmetric theories automatically include gravity [^]. Un- 
fortunately supergravity models are inherently non-renormalizible, which 
prevents up to now clear predictions. Nevertheless, the spontaneous sym- 
metry breaking of supergravity is important for the low energy spectrum 
of supersymmetry [p9|] . The most common scenario is the hidden sector 



scenario 



40|] , in which one postulates two sectors of fields: the visible sector 
containing all the particles of the GUT's described before and the hidden 
sector, which contains fields which lead to symmetry breaking of supersym- 
metry at some large scale Asusy- One assumes that none of the fields in the 
hidden sector contains quantum numbers of the visible sector, so the two 
sectors only communicate via gravitational interactions. Consequently, the 
effective scale of supersymmetry breaking in the visible sector is suppressed 
by a power of the Planck scale, i.e. 

Msusy - (4.6) 

where n is model-dependent (e.g. n = 2 in the Polonyi model). Thus the 



37 



Figure 4.3: Examples of proton decay in the minimal supersymmetric model via 
wino and Higgsino exchange. 

SUSY breaking scale can be large, above 10^° GeV, while still producing 
a small breaking scale in the visible sector. In this case the fine-tuning 
problem can be avoided in a natural way and it is gratifying to see that the 
first experimental hints for MgusY are indeed in the mass range consistent 
with eq. 



The hidden sector scenario leads to an effective low-energy theory with 
explicit soft breaking terms, where soft implies that no new quadratic di- 



vergences are generated [^. The soft-breaking terms in string- inspired su- 
pergravity models have been studied recently in refs. [Q. A final theory, 
which simultaneously solves the cosmological constant problem and ex- 



plains the origin of supersymmetry breaking, needs certainly a better un- 
derstanding of superstring theory. 

The unification scale in SUSY is large 



As discussed in chapter the limits on the proton lifetime require the 
unification scale to be above 10^^ GeV, which is the case for the MSSM. 
In addition, one has to consider proton decay via graphs of the type shown 



in fig. ^4.3| . These yield a strong constraint on the mixing in the Higgs 



sector as will be discussed in detail in the last chapter. 



Prediction of dark matter 

The lightest supersymmetric particle (LSP) cannot decay into normal mat- 
ter, because of R-parity conservation (see the next section for a definition 
of R-parity). In addition R-parity forbids a coupling between the LSP and 
normal matter. 



Consequently, the LSP is an ideal candidate for dark matter [^ , which is 
believed to account for a large fraction of all mass in the universe (see next 
chapter). The mass of the dark matter particles is expected to be below 
one TeVi 



4.2 SUSY interactions 

The quantum numbers and the gauge couplings of the particles and sparticles 
have to be the same, since they belong to the same multiplet structure. 

The interaction of the sparticles with normal matter is governed by a new mul- 
tiplicative quantum number called R-parity, which is needed in order to prevent 



38 



baryon- and lepton number violation. Remember that quarks, leptons and Hig- 
gses are all contained in the same chiral supermultiplet, which allows couplings 
between quarks and leptons. Such transitions, which could lead to rapid proton 
decay, are not observed in nature. Therefore, the SM particles are assigned a pos- 
itive R-parity and the supersymmetric partners are R-odd. Requiring R-parity 
conservation implies that: 

• sparticles can be produced only in pairs 

• the lightest supersymmetric particle is stable, since its decay into normal 
matter would change R-parity. 

• the interactions of particles and sparticles can be different. For example, 
the photon couples to electron-positron pairs, but the photino does not 
couple to selectron-spositron pairs, since in the latter case the R-parity 
would change from -1 to +1. 



4.3 The SUSY Mass Spectrum 

Obviously SUSY cannot be an exact symmetry of nature; or else the supersym- 
metric partners would have the same mass as the normal particles. As mentioned 
above, the supersymmetric partners should be not too heavy, since otherwise the 
hierarchy problem reappears. 

Furthermore, if one requires that the breaking terms do not introduce quadratic 
divergences, only the so-called soft breaking terms are allowed [p] ]. 

Using the supergravity inspired breaking terms, which assume a common mass 
mi/2 for the gauginos and another common mass mo for the scalars, leads to the 
following breaking term in the Lagrangian (in the notation of ref. 



^Breaking = — "^Q X! I I ^ ~ "^1/2 X! "^a-^a 
i a 

- Amo [K,QaU',H2 + K,Q,DIH^ + K^L^EIH,] - Bm^ [fiH.m]^) 

Here 

h^f''^ are the Yukawa couplings, a,b = 1,2,3 run over the generations 
Qa are the SU(2) doublet quark fields 

are the SU(2) singlet charge-conjugated up-quark fields 
are the SU(2) singlet charge-conjugated down-quark fields 
La are the SU(2) doublet lepton fields 

are the SU(2) singlet charge-conjugated lepton fields 
Hi^2 are the SU(2) doublet Higgs fields 
(pi are all scalar fields 
\a are the gaugino fields 
The last two terms in C Breaking originate from the cubic and quadratic terms 
in the superpotential with A, B and /x as free parameters. In total we now have 
three couplings ctj and five mass parameters 

mo, mi/2, nit), A{t), Bit) 
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with the following boundary conditions at Mqut (t = 0): 

scalars : mg = = m|, = m| = m|; = ml; (4-9) 
gauginos : Mi = mi/2, i = 1,2,3; (4-10) 
couplings : «j(0) = acuT, i = 1,2,3. (4-11) 

Here Mi, M2, and M3 are the gauginos masses of the U{1), SU{2) and SU{3) 
groups. In = 1 supergravity one expects at the Planck scale B = A — 1. With 
these parameters and the initial conditions at the GUT scale the masses of all 
SUSY particles can be calculated via the renormalization group equations. 



4.4 Squarks and Sleptons 

The squark and slepton masses all have the same value at the GUT scale. How- 
ever, in contrast to the leptons, the squarks get additional radiative corrections 
from virtual gluons (like the ones in fig. |3.2| for quarks), which makes them heavier 
than the sleptons at low energies. These radiative corrections can be calculated 
from the corresponding RGE, which have been assembled in the appendix. The 
solutions are: 



{t = 


66) 


= ml 


+ 0.52m?/2 - 0.27cos(2/3)M| 


(4.12) 


KSt = 


66) 


= ml 


+ 0.52m2/2 + 0.5 cos(2/?)M| 


(4.13) 


ml^{t = 


66) 


= ^1 


+ 0.15m?/2 - 0.23cos(2/5)M| 


(4.14) 


ml^{t = 


66) 


= ml 


+ 6.6m2/2 + 0.35 cos(2/?)M| 


(4.15) 


ml^{t = 


66) 


= ml 


+ 6.6m2/2 - 0.42 cos(2/5)M| 


(4.16) 




66) 


= ml 


+ 6.2m?/2 + 0.15 cos(2/3)M| 


(4.17) 


ml^{t = 


66) 


= ml 


+ 6.1m?/2 - 0.07cos(2/5)M|, 


(4.18) 



where /3 is the mixing angle between the two Higgs doublets, which will be de- 
fined more precisely in section O. The coefficients depend on the couplings as 
shown explicitly in the appendix. They were calculated for the parameters from 
the typical fit shown in table §]T] {acuT = 1/24.3, Mqut = 2.0 ■ 10^*^ GeV and 
sin^ 9w = 0.2324). For the third generation the Yukawa coupling is not necessar- 
ily negligible. If one includes only the correction from the top Yukawa coupling 
YtQ, one finds: 

(4.19) 
(4.20) 
(4.21) 
(4.22) 



mlit = 


66) 


= 






mUt = 


66) 


= 


- 0.4< 


— 1.2\m\i2 


Kit = 


66) 


= ^Ir 


+ "^? 


- 0.96mo - 2.42mi/2 


K(t = 


66) 


= ^II 




- 0.48mo - 1.21mi/2 



-'^For large values of the mixing angle tan (3 in the Higgs sector, the b-quark Yukawa coupling 
can become large too. However, since the limits on the proton lifetime limit tan (3 to rather 
small values (see last chapter), this option is not further considered here. 
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The numerical factors have been calculated for At{0) = and the explicit depen- 
dence on the couplings can be found in the appendix. Note that only the left- 
handed b-quark gets corrections from the top-quark Yukawa coupling through a 
loop with a charged Higgsino and a top-quark. The subscripts L or i? do not 
indicate the helicity, since the squarks and sleptons have no spin. The labels 
just indicate in analogy to the non-SUSY particles, if they are SU{2) doublets 
or singlets. The mass eigenstates are mixtures of the L and R weak interaction 
states. Since the mixing is proportional to the Yukawa coupling, we will only 
consider the mixing for the top quarks. After mixing the mass eigenstates are: 
(using the same numerical input as for the light quarks): 



Ut = 66) 



1 

2 
1 

2 
1 

±- 
2 



0.6ml + 9.2mi/2 + 2m'^t - 0.19 cos(2/3)M| 



1.6m^ 



/2 



0.5ml - 0-5 cos(2/3)M|] + 4m|( + /// tan/3)2 

(4.23) 



where the values of At and /i at the weak scale can be calculated as: 



At{Mz) = 4.6A(0) + 1.7 
^(Mz) = 0.63/i(0) 



mi/2 
mo 



(4.24) 
(4.25) 



Note that for large values of At{0) or fi combined with a small tan/3 the 
splitting becomes large and one of the stop masses can become very small; since 



the stop mass lower limit is about 46 GeV[0, this yields a constraint on the 
possible values of mt, mi/2, tan/5 and /x. 



4.5 Charginos and Neutralinos 

The solutions of the RGE group equations for the gaugino masses are simple: 

M^{t) = ?^mv2. (4.26) 



Numerically at the weak scale (t = 2\n{MGUT/Mz) = 66) one finds (see fig. [4.4|) : 



Msig) ^ 2.7mi/2, (4.27) 
M2iMz) ^ 0.8mi/2, (4.28) 
Mi(Mz) ^ 0.4mi/2. (4.29) 

Since the gluinos obtain corrections from the strong coupling constant a^, they 
grow heavier than the gauginos of the SU{3)c ® SU{2)l ® U{1)y group. 

The calculation of the mass eigenstates is more complicated, since both Hig- 
gsinos and gauginos are spin 1/2 particles, so the mass eigenstates are in general 
mixtures of the weak interaction eigenstates. The mixing of the Higgsinos and 
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Figure 4.4: Typical running of the squark (q), slepton (bl), and gaugino 
(Ml, M2, Ms) masses (solid lines). The dashed hues indicate the running of 
the four neutralinos and two charginos. 
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gauginos, whose mass eigenstates are called charginos and neutralinos for the 
charged and neutral fields, respectively, can be parametrised by the following 
Lagrangian: 



^Gaugino—Higgsino 

where Aq, a = 1, 2, . . . , 8, are the Majorana gluino fields and 



X 



( B \ 

Hi 




are the Majorana neutralino and Dirac chargino fields, respectively. Here all the 
terms in the Lagrangian were assembled into matrix notation (similarly to the 
mass matrix for the mixing between B and in the SM, eq. pT7|) . The mass 



matrices can be written as 14 : 



Ml 
Ms 

Mz cos /? sin Q-^ Mz cos fi cos Qw 

V Mz sin /? sin dw —Mz sin /? cos Qw 



Ma 

^M^Jv cos /5 



-Mz cos /3 sin Q-w Mz sin /5 sin Q-w \ 
Mz cos /3 cos 6'vi/ — M^ sin /? cos 6*;^ 
-/i 
-/i 

(4.30) 

v^MH/sin/5^ ^^g^^ 



The last matrix leads to two chargino eigenstates Xi 2 with mass eigenvalues 



Ml. 



M^ + + 2M^ T V (^2 -/"=')' + 4M^ cos2 2/3 + 4M2,(M| + /i^ + 2M2/i sin 2/3) 

(4.32) 



The dependence on the parameters at the GUT scale can be estimated by sub- 
stituting for M2 and /x their values at the weak scale: M2(M^) ^ 0.8mi/2 and 
li{Mz) ~ 0.63/i(0). In the case favoured by the fit discussed in chapter || one 
finds ^ » M2 ~ Mz, in which case the charginos eigenstates are approximately 
M2 and fi. 

The four neutralino mass eigenstates are denoted by x°(i = 1,2,3,4) with 
masses M^o < ■ ■ ■ < M^o. The sign of the mass eigenvalue corresponds to the 
CP quantum number of the Majorana neutralino state. 

In the limiting case Mi,M2,/i >> Mz one can neglect the off-diagonal ele- 
ments and the mass eigenstates become: 



Xi 



[5,1^3, ^(^1-^2), ^(^1 + ^2)] 



(4.33) 



with eigenvalues |Mi|, IM2I, and respectively. In other words, the bino 
and neutral wino do not mix with each other nor with the Higgsino eigenstates 
in this limiting case. As we will see in a quantitative analysis, the data indeed 
prefer Mi, M2, /x > Mz, so the LSP is bino- like, which has consequences for dark 
matter searches. 
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4.6 Higgs Sector 

The Higgs sector of the SUSY model has to be extended with respect to the one 
of the SM for two reasons: 

• the Higgsinos have spin 1/2, which imphes they contribute to the gauge 
anomaly, unless one has pairs of Higgsinos with opposite hypercharge, so 
in addition to the Higgs doublet with Yw=l one needs a second one with 
Yw=-l: 



^i(l,2,-l)= (fi j, i/2(l,2,l)=(^i j (4.34) 

• The introduction of the second Higgs doublet solves simultaneously the 
problem that a single doublet can give mass to only either the up- or down- 
type quarks, as is apparent from the fact that only the neutral components 
have a non zero vev, since else the vacuum would not be neutral. So one 
can write: 

<^.>=(';). <^^>=(°)- ("5) 

In the SM the conjugate field can give mass to the other type. However, 
supersymmetry is a spin-symmetry, in which the matter - and Higgs fields 
are contained in the same chiral supermultiplet. This forbids couplings 
between matter fields and conjugate Higgs fields. With the two Higgs fields 
introduced above. Hi generates mass to the down-type matter fields, while 
H2 generates mass for the up-type matter fields. 

The supersymmetric model with two Higgs doublets is called the Minimal Su- 
persymmetric Standard Model (MSSM). The mass spectrum can be analyzed by 
considering again the expansion around the vacuum expectation value, given by 



eq. 4.35 



( vi + ^ (H^ cos a-h°sma + sin p - iG^ sin p) \ . . 

^ H-sin(3-G'cosf3 ) ^^'^'^ 

\ V2 + ^ (H^ sin a + h'^ cos a + iA'^ cos p + iG'^ sin P) J 

Here H, h and A represent the fluctuations around the vacuum corresponding to 
the real Higgs fields, while the G's represent the Goldstone fields, which disappear 
in exchange for the longitudinal polarization components of the heavy gauge 
bosons. The imaginary and real sectors do not mix, since they have different 
CP-eigenvalues; a and P are the mixing angles in these different sectors. The 
mass eigenvalues of the imaginary components are CP-odd, so one is left with 
2 neutral CP-even Higgs bosons and h^, 1 CP-odd neutral Higgs bosons 
and 2 CP-even charged Higgs bosons. 
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The complete tree level potential for the neutral Higgs sector, assuming colour 
and charge conservation, reads: 



(4.38) 

Note that in comparison with the potential in the SM, the first terms do not have 
arbitrary coefficients anymore, but these are restricted to be the gauge coupling 
constants in supersymmetry, again because the Higgses belong to the same chiral 
multiplet as the matter field^. The last three terms in the potential arise from 
the soft breaking terms with the following boundary conditions at the GUT scale: 

ml{0) = ml{0) = fi{OY + ml ml{0) = -Bfi{0)mo, (4.39) 

where /i(0) is the value of /i at the GUT scale. Since /i generates mass for the 
Higgsinos, one expects /i to be small compared with the GUT scale. A low /i value 
can be obtained dynamically, if one adds a singlet scalar field to the MSSM. This 
will not be considered further. Instead fi is considered to be a free parameter to 
be determined from data (see chapter H). 

From the potential one can derive easily the five Higgs masses in terms of 
these parameters by diagonalization of the mass matrices: 

2 ^ 1 9'Vh . . 

2d(f),d(Pj ^ ' 

where (^j is a generic notation for the real or imaginary part of the Higgs field. 
Since the Higgs particles are quantum field oscillations around the minimum, 
eq. |4.4CI| has to be evaluated at the minimum. One finds zero masses for the 



Goldstone bosons. These would-be Goldstone bosons and are "eaten" by 
the SU{2) gauge bosons. For the masses of the five remaining Higgs particles one 
finds Hq: 

i neutral Higgs A: 

m\ = ml + ml. (4.41) 
ml± = m\ + M^. (4.42) 



CP-odd neutral Higgs A: 
Charged Higgses H^: 
CP-even neutral Higgses if, h: 



1 



H.h 



m\ + Ml± J{m\ + M|)2 - 4m^M| cos^ 2/3 



(4.43) 



2 

By convention itih > iTih- The mixing angles a and /3 are related by 

fjl'^ _)_ j^^2 

tan 2a = ^ f tan 2/3 (4.44) 



^In principle one should consider the running of the gauge couplings between the electroweak 
scale and the mass scale of the Higgs bosons. However, since the Higgs bosons are expected to 
be below the TeV mass scale, this running is small and can be neglected, if one considers all 
other sources of uncertainty in the MSSM. 
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Figure 4.5: Corrections to the Higgs self-energy from Yukawa type interactions. 

Vi and V2 have been chosen real and positive, which implies < (3 < 7t/2. 
Furthermore, the electroweak breaking conditions require tan/5 > 1, so 

7r/4 < /3 < 7i/2. (4.45) 

From the mass formulae at tree level one obtains the once celebrated SUSY mass 
relations: 

mH±>Mw (4.46) 

mh<mA< Mh (4.47) 

rrih < Mz cos 2f3 < Mz (4.48) 

ml + m]j = m\ + Ml. (4.49) 

After including radiative corrections the lightest neutral Higgs mh becomes con- 
siderably heavier and these relations are not valid anymore. The mass formulae 
including the radiative corrections are given in the appendix. 



4.7 Electroweak Symmetry Breaking 

The coupling yU plays an important role in the shape of the potential and con- 
sequently in the pattern of electroweak symmetry breaking, which occurs if the 
minimum of the potential is not obtained for < Hi >=< H2 >=0. In the 
SM this condition could be introduced ad-hoc by requiring the coefficient of the 
quadratic term to be negative. In supersymmetry this term is restricted by the 



gauge couplings ||34||. A non-trivial minimum can only be obtained by the soft 
breaking terms, if the mass matrix for the Higgs sector, given by Mfj = q^.qjj. , 
has a negative eigenvalue. This is obtained if the determinant is negative, i.e. 

\ml{t)\'^ > ml{t) ml{t). (4.50) 

In order that the new minimum is below the trivial minimum with< Hi >=< 
H2 >=0, one has to require in addition Vh{vi,V2) < Vh{0, 0) < Vh{oo, 00), which 
is fulfilled if 

ml{t) + mlit) > 2\ml{t)\. (4.51) 
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If one compares eqns. [4.50| and [4.51j and notices from eq. [4.39| that mi = m2 



one realizes that these conditions cannot be fulfilled simultaneously, at least not 
at the GUT scale. 

However, at lower energies there are substantial radiative corrections, which 
can cause differences between m2 and mi, since the first one involves mass cor- 
rections proportional to the top Yukawa coupling Yt{0), while for the latter these 
corrections are proportional to the bottom Yukawa coupling. Typical diagrams 
are shown in fig. |4.5| . From the RGE for the mass parameters in the Higgs 
potential one finds at the weak scale: 

^^{t = 66) = 0.40/1^(0) (4.52) 
ml{t = 66) = m^ + 0.40/2^(0) + 0.52m2/2 (4.53) 

ml{t = 66) = -0.44m2 + 0.40/^2(0) - 3. llm2/2 

-0.09At(0)momi/2 - 0.02A(0)^mg (4.54) 
ml{t = 66) = 0.63m^(0) + 0.04/i(0)mi/2 + 0.19A(0)moAi(0). (4.55) 

The coefficients were evaluated for the parameters of the fit to the experimental 
data (central column of table |6]^ in chapter The explicit dependence of the 



coefficients on the coupling constants is given in the appendix. The coefficients of 
the last three terms in m2 depend on the top Yukawa coupling. This dependence 
disappears if the masses of the stop and top quarks in the diagrams of fig. 
are equal. However, if the stop mass is heavier, the negative contribution of the 
diagram with the top quarks dominates; in this case m2 decreases much faster 
than mi with decreasing energy and the potential takes the form of a mexican 
hat, as soon as conditions [4.50| and [4.51| are satisfied. Since A{t) is expected 



to be small, the dominant negative contribution is proportional to mi/2 (see eq. 



4.54| ), so the electroweak breaking scale is a sensitive function of both the initial 
conditions, the top Yukawa coupling and the gaugino masses. 
The minimum of the potential can be found by requiring: 

9V 2 o 2 , 9^ + 9"^f 2 n 
2m-|^fi — 2m^V2 H ("^i — V2)Vi = (J 



- 2mlv2-2mlv,-l±^{vl-vl)v2 = ^ 



Here we substituted 

< Hi >= Vl = V cos P, < H2 >= V2 = V sin /3, 

where 

From the minimization conditions given above one can derive easily: 

.2 ^ 1 

(/ + ^?'2)(tan^/3 
2ml = (ml + ml)sm2P (4.57) 
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Ml^l±^v- = (4.58) 

M^ = |-^;2 = Mlcos'^Ow (4.59) 

The derivation of these formulae including the one-loop radiative corrections 
is given in the appendix. 
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Chapter 5 

The Big Bang Theory 



5.1 Introduction 

In the 1920's Hubble discovered that most galaxies showed a redshift in the visible 
spectra, implying that they were moving away from each other. This observation 
is one of the basic building blocks of the Big Bang Theory [Q, which assumes the 
universe is expanding, thus solving the problem that a static universe cannot be 
stable according the Einstein's equations of general relativity. An expanding 
universe will cool down, so at the beginning the universe might have been hot. 
The remnants of the radiation of such a hot universe can still be observed today 
as microwave background radiation corresponding to a temperature of a few de- 
grees. This radiation was first predicted by Gamow, but accidentally observed 
in 1963 by Penzia and Wilson from Bell Laboratories]^ as noise in microwave 
antennas used for communication with early satellites. Such an antenna is only 
sensitive to a single frequency. Recently, the whole spectrum was measured by 
the COBE0 satellite and it was found to be indeed describable by a black body 
radiation, as shown in fig. |5.1| (from ref. 0]). The deviation from perfect isotropy. 



if one ignores the dipole anisotropy from the Doppler shift caused by the move- 
ment of the earth through the microwave background, is a few times 10~^. This 
has strong implications for theories concerning the clustering of galaxies, since 
this radiation was released soon after the "bang" and hardly interacted after- 
wards, so the inhomogeneities in this radiation are proportional to the density 
fluctuations in the early universe! These density fluctuations are the seeds for 
the final formation of galaxies. As will be discussed later, such small anisotropics 
have strong implications for the models trying to understand the formation of 
galaxies and the nature of the dark matter in the universe. Direct evidence that 
the temperature from the microwave background is indeed the temperature of 
the universe came from the measurement of the temperature of gas clouds deep 
in space. As it happens, the rotational energy levels of cyanogen {CN) are such 
that the 3K background radiation can excite these molecules. From the detection 
of the relative population of the groundstate and the higher levels the excitation 



temperature was determined to be Tcn = '^-^'^^-omi K||5^, which is in excel 



^They were awarded the Nobel prize for this discovery in 1978. 
^Cosmic Background Explorer. 
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Frequency (cycles/cm) 



Figure 5.1: Spectrum of the microwave background radiation as measured by 
the COBE satelhte. The curve is the black body radiation corresponding to a 
temperature of 2.726 K. 



lent agreement with the direct measurement of the microwave background of 
2.726 ± 0.010 K by COBEpTl. 



Other evidence that the universe was indeed very hot at the beginning came 
from the measurement of the natural abundance of the light elements: the uni- 
verse consists for 74% out of hydrogen, 24% helium and 1% for the remaining 
elements. Both the rarity of heavy elements and the large abundance of helium 
are hard to explain, unless one assumes a hot universe at the beginning. 

The reason for the low abundance of the heavier elements in a hot universe 
is simple: they are cracked by the intense radiation around. The abundance of 
the light elements is plotted in fig. |5.2| as function of the ratio rj of primordial 
baryons and photons (from ref. [0). Agreement with experimental observations 
can only be obtained for t] in the range 3 — 7 ■ 10^^°. 

The very heavy elements can be produced only at much lower temperatures, 
but high pressure, so it is usually assumed that the heavy elements on earth and 
in our bodies were cooked by the high pressure inside the cores of collapsing stars, 
which exploded as supernovae and put large quantities of these elements into the 
heavens. They clustered into galaxies under the influence of gravity. 

The ratio of helium and hydrogen is determined by the number of neutrons 
available for fusion into deuterium and subsequently into helium at the freeze-out 
temperature of about 1 MeV or 10^° K.At these high temperatures no complex 
nuclei can exist, only free protons and neutrons. They can be converted into 
each other via charged weak interactions like ep ^ Ugn and en ^ V^p. Note 
that this is the same interaction which is responsible for the decay of a free 
neutron into a proton, electron and antineutrino. The weak interactions maintain 



50 



T I I M I I 




Figure 5.2: Big Bang nucleosynthesis predictions for the primordial abundance of 
the light elements as function of the primordial ratio rj of baryons and photons. 
From ref. [49]. 



thermal equilibrium between the protons and neutrons as long as the density and 
temperature are high enough, thus leading to a Boltzmann distribution: 

- = e-^/'^^, (5.1) 
P 

where Q = (m„ — mp)c^ = 1.29 MeV is the energy difference between the states. 
Thermal equilibrium is not guaranteed anymore if the weak interaction rates F 
are slower than the expansion rate of the universe given by the Hubble constant, 
i.e. freeze-out occurs when F < H(t). This happens by the time the temperature 
is about 10^ K or 0.1 MeV. Then the ratio n/p is about 1/7. Since the photon 
energies at these temperatures are too low to crack the heavier nuclei, nuclei can 
form through reactions hke n+p ^H+p and '^H+n 

which in turn react to form ^i7e. The latter is a very stable nucleus, which hardly 
can be cracked, so the chain essentially stops till all neutrons are bound inside 
^iJe! Heavier nuclei are hardly produced at this stage, since there are no stable 
elements with 5 or 8 nuclei, so as soon as ^He catches another nuclei it will decay. 
Consequently the n/p ratio 1/7, as determined from the Boltzmann distribution, 
yields a '^He mass fraction Y^g 

= 4^ « I- (5.2) 

n/p + 1 4 
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Figure 5.3: The lineshape for different number of neutrino types. The data 
(black points) exclude more than three types with mass below Mz/2 ^ 45 GeV. 

Experimentally the mass fraction Yp of ^ife is 23 ± 1% |5^! The mass fractions 
of deuterium, ^He, and '^Li are many orders of magnitude smaller|[l|, The 
concentration of the latter elements is a strong function of the primordial baryon 
density (see fig. |5.2| ), since at high enough density all the deuterium will fuse into 
'^He, thus eliminating the "components" for ^He and '^Li. 

As said above, freeze-out occurs, if F < H{t). Thus the expansion rate H{t) 
around T ^ 1 MeV determines the ^ife abundance. The expansion rate in turn is 
determined by the fraction of relativistic particles, like neutrinos, light photinos 
etc. Roughly for each additional species the primordial ^ife abundance increases 



by 1%, as shown in fig. 5^ for a neutron half-life time of 10.6 minutes (from ref. 
^)|. The neutron lifetime is not negligible on the scale of the first three minutes, 
so it has to be taken into account. If one wants to reconcile the abundance of 
all light elements, there can only be three neutrino species (see fig. |5.2| ) with 
practically no room for other weakly interacting relativistic particles like light 
photinos! Present collider data confirm that there are indeed only three light 
neutrinos [pO| 



= 2.99 ± 0.04. (5.3) 

The strongest constraint comes from the resonance data at LEP[Q. An 
example of the quality of the data is shown in fig. p.3| . Note that collider data 
limit the number of neutrino generations, while nucleosynthesis is sensitive to 
all kinds of light particles, where light means about one MeV or less. Happily 

■^The presently accepted value is 10.27 ± 0.024 minutes[B3l. 
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enough the colhder data require the hghtest neutrahno to be above 18.4 GeVpO 
so there is no conflict between cosmology and supersymmetry. 

Note that this is a beautiful example of the interactions between cosmology 
and elementary particle physics: from the Big Bang Theory the number of rela- 
tivistic neutrinos is restricted to three (or four if one takes the more conservative 
upper limit on the ^He abundance to be 0.25) and at the LEP accelerator one 
observes that the number of light neutrinos is indeed three! Alternatively, one 
can combine the accelerator data and the abundance of the light elements to 
"postdict" the primordial helium abundance to be 24% and use it to obtain an 
upper limit on the baryonic density [^ : 

Pb < O.lpc, (5.4) 

where pc is the critical density needed for a flat universe. The critical density 
will be calculated in section p.4| . Thus LEP data in combination with baryoge- 
nesis strengthens the argument that we need non-baryonic dark matter in a fiat 
universe, for which p = p^. Other arguments for dark matter will be discussed in 
section |5.12 . 



In spite of the marvelous successes of this model of the universe, many ques- 
tions and problems remain, as mentioned in the Introduction. However, GUT's 
can provide amazingly simple solutions, at least in principle, since many details 
are still open. 

These problems will be discussed more quantitatively in the next sections, 
starting with Einsteins equations in a homogeneous and isotropic universe, con- 
ditions which have been well verified in the present universe and which make the 
solutions to Einstein's equations particularly simple. Especially, it is easy to see 
that a phase transition can lead to inflation, the key in all present cosmological 
theories. 



5.2 Predictions from General Relativity 

At large distances the universe is homogeneous, i.e. one finds the same mass 
density everywhere in the universe, typically 

Puniv = (4 - 16) 10-27 kg/m^ (5.5) 

which corresponds to 2.5 -10 hydrogen atoms per cubic meter. (In comparison, an 
extremely good vacuum of lQ~'^N/m'^ at 300 K contains about 2 10^^ molecules 
per cubic meter. Of course, the volume to be averaged over should be chosen to be 
much larger than the size of clusters of galaxies. Furthermore, the same density 
and temperature is observed in all directions, i.e. the universe is very isotropic. If 
no point and no direction is preferred in the universe, the possible geometry of the 
universe becomes very simple: the curvature has to be the same everywhere, i.e. 
instead of a curvature tensor one needs only a single number, usually written as 
K{t) = k/ R^{t), where R{t) is the so-called scale factor. This factor can be used 
to define dimensionless time-independent (comoving) coordinates in an expanding 
universe: the proper (or real) distance D{t) between two galaxies scales as 

D{t) = R{t)d, (5.6) 
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where d is the distance at a given time to- The factor k introduced above defines 
the sign of the curvature: k = imphes no curvature, i.e. a fiat universe, while 
k = 1) corresponds to a space with a positive curvature (spherical) and 

k = —1 corresponds to a space with a negative curvature (hyperbolic). 

The movement of a galaxy in a homogeneous universe can be compared to the 
molecules in a gas; the stars are just the atoms of a molecule and the molecules 
are homogeneously distributed. Differentiating equation EM results in 



V = R{t)d, (5.7) 
or substituting d from eq. |5.(j| results in the famous relation between the velocity 



and the distance of two galaxies: 

m 

R{t) 



D{t) = H{t)D{t), (5.8) 



where H{t) is the famous Hubble constant. 

This relation between the velocity and the distance of the galaxies was first 
observed experimentally by Hubble in the 1920's. He observed that all neigh- 
bouring galaxies showed a redshift in the spectral lines of the light emitted by 
specific elements and the redshift was roughly proportional to the distance. So 
this was the first evidence that we are living in an expanding universe, which 
might have been created by a "Big Bang" . 

The Hubble relation |5.8| is a direct consequence of the homogeneity and 
isotropy of the universe, since the scale factor cannot be a constant in that case. 
This follows directly from Einstein's field equations of general relativity, which 
can be written as: 

Rit) = -^Ht) + Mt))R{t), (5.9) 



R{ty 8n G , , kc^ 

7#)-ljr"W = -ip. (5.10) 

where G = 6.67 • 10~^^Nm^/kg^ is the gravitational constant, R and R are the 
derivatives of R with respect to time, p is the pressure and u is the energy density. 

A static universe, in which the derivatives and pressure are zero, implies u{t) = 
0, so a static universe cannot exist unless one introduces additional potential 
energy in the universe, e.g. Einstein's cosmological constant. At present there is 



no experimental evidence for such a term[^ . 



5.3 Interpretation in terms of Newtonian Me- 
chanics 

Since the energy density, and correspondingly the curvature, is small in our 



present universe, relativistic effects can be neglected and the field equations 5.£ 



and 5.10 have a simple interpretation in terms of Newtonian mechanics. Consider 
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a spherical shell with radius R and mass m. The mass inside this sphere can be 
related the average density p: 



3 ^ 



(5.11) 



For an expanding universe the total mechanical energy of the mass shell can be 
written as the sum of the kinetic and potential energy: 



tot 



1 -9 GMm 
-mW 

2 R 



-mR' 
2 



R^ 8 

— vrGp 

i?2 3 



(5.12) 



The expression in brackets is just the left hand side of eq. ^.10| , so the sum of 
kinetic - and potential energy determines the sign of the curvature k. \{ k = 1, 
then Etot < 0, implying that the universe will recollapse under the influence of 
gravity ("Big Crunch"), just like a rocket which is launched with a speed below 
the escape velocity, will return to the earth; /c = — 1, on the other hand, implies 
that the universe will expand and cool forever ( "Big Chill" ) . In case of /c = the 
total energy equals zero (flat Euclidean space), in which case the gravitational 
energy, or equivalently the mass density, is sufficient to halt the expansion. 

The first field equation (eq. |0| ) follows from the second equation (eq. |5.10|) 
by differentiation and taking into account that in an expanding universe energy 
is converted into gravitational potential energy: when the volume increases by 
an infinitesimal amount AV^, then the remaining energy in the gas decreases by 
an amount p/SVphys^ where p denotes the pressure. Therefore 



P(t)Vphys{t). 



(5.13) 



Here we used Vphys = Vo-R^(t), where Vq is the volume in comoving time- independent 



coordinates analogous to eq. |5.6| . On the other hand follows from E{t) 

u{t)Vphys{t) 



E{t) = u{t)Vpr,y,{t) + ^^^u{t)Vp,,ys{t). 



Combining eqs. |5.13| and |5.14| results in: 



u{t) 



u{t)+p{t)). 



(5.14) 



(5.15) 



Substituting this equation for u{t) after differentiation of eq. |5.10| yields eq. O 



5.4 Time Evolution of the Universe 

To find out how the universe will evolve in time, one needs to know the equa- 
tion of state, which relates the energy density to pressure. Usually energy and 
pressure are proportional, i.e. p = apc^, where a = for cold non-relativistic 
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Figure 5.4: Evolution of the radius of the universe for a closed (C), flat (F) or 
open (O) universe with and without inflation. From ref. [3]. 



matter {p = 0) and a = 1/3 for a relativistic hot gas, as follows from elementary 
Thermodynamics. From eq. |5.15| given above, it follows immediately that 



P 



(5.16) 



and substituting this into eq. 5.10 results in: 



R oc t3(i+"). 



(5.17) 



if we neglect the curvature term, i.e. either R is large or k small. As we will see, 
both are true after the inflationary phase of the universe. 



hot relativistic 


t > tl 


R oc t^' 


P = Di/R' 


inflation 


tl > t> t2 


i?oc e^* 


p = const. 


hot relativistic 


t2>t>ts 


R oc 


P = D2/R^ 


cold non — relativ. 


t>t3 


R oc ^2/3 


P = D,/R^ 



Table 5.1: The time dependence of the scale factor and energy density during 
various stages in the evolution of the universe. Typically, ti ~ lO"'^^ s, ^2 ~ 10~^^ 
s, and ts ~ 10^ yrs. The constants Di are integration constants. 



The time dependence has been summarized in table ^.1| for various stages of 
the universe. The inflationary period will be discussed in the next section. One 
observes that the scale factor vanishes at some time t = and the energy density 
becomes infinite at that time. This singularity explains the popular name "Big 
Bang" theory for the evolution of the universe. The solutions for R{t) are shown 



graphically in fig. |5.4| for three flat universe (fc = 0), an open universe 



{k = —1) and a closed universe {k = 1) (from ref.| 
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An open or flat universe will expand forever, since the kinetic energy is larger 
than the gravitational attraction. A closed universe will recollapse. The lifetime 
of a closed universe with p > —p/3 and a total mass M of cold non-relativistic 
matter is0: 

3 Mp ' ^ ^ 

so the present lifetime of the universe of at least 10^° yrs gives a strong upper 
limit on the density of the universe. 

The lifetime can easily be calculated, if we assume a flat universe: from table 
O it follows that R{t) oc t^/^ for most of the time. Substituting this and its time 



derivative into the deflnition of the Hubble constant (eq. |5]^ ) results in: 

H{t) = ^. (5.19) 
With the presently accepted value of the measured Hubble constant: 

fcTD 

H = 100 ho (— — ) ^ K (3 10'')-h-' ^ ho 10-1° yrs-\ (5.20) 
s Mpc 

where h indicates the experimental uncertainty (0.4 < ho < 1), one flnds for the 
age of the universe: 

tuniverse = 2/3H = 2/(3/io) " lO^Vrs. (5.21) 

The critical density, which is the density corresponding to a flat universe, can 
be calculated from eqn. ^.12| by requiring Etot = (or equivalently k = 0)and 
substituting for R/R the Hubble constant (see eq. 



3^2 

= ^ = 2 ■ 10"'' hi kg/m\ (5.22) 



where the numerical value of H from eq. |5.20| was used. 

The size of the observable universe, the horizon distance D^, can be calculated 
in the following way: the proper distance between two points is R(t)d (eq. p.6| ), 
where d is the distance in comoving coordinates. Light propagates on the light- 
cone. This can be studied most easily by considering the time t] in comoving 
coordinates with 

dt = R{t)dr]. (5.23) 

In comoving coordinates the distance dh light can propagate is cdrj, so dh = 
cJdT] = cJ dt/R{t) or the proper distance Dh = R{t)dh equals: 

For most of the time R{t) = at^^^ (see table |0| ). Substituting this into eq. 
yields: 

Dh = 3ct = 2c/H{t) = 0.9ho^ ■ 10^^ m, (5.25) 



where for the lifetime t of the universe eqs. 5. IS and 5.20 were used. 
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5.5 Temperature Evolution of the Universe 



In the previous section the scale factor and the energy density were calculated 
as function of time. The energy density has two components: the energy density 
from the photon radiation in the microwave background prad and the energy 
density of the non-relativistic matter Pmatter- At present p^ad is negligible, but 
at the beginning of the universe it was the dominating energy. Assuming this 
radiation to be in thermal equilibrium with matter implies a black body radiation 
with a frequency distribution given by Planck's law and an energy density 

prad = aT^, (5.26) 



where a = 7.57 ■ 10 J m ^. Since prad oc l/i?^(see table one finds from 
eq. p. 26: 



T oc l/i?(t), (5.27) 

from which follows immediately: R/R = —T/T and R~^ oc T^. Substituting 
these expressions and eq. p.26| into the second field equation (eq. p.lO| ) leads to: 

-T^T\ 5.28 

since the term with kc^ is only proportional to T^, so it can be neglected at high 
temperatures. Integrating eq. |5.28| yields: 

1/4 



3c^ \ ' 1 . . . .0 /Is _ /I s 



^ = 1.5 • 10^'^ K-d — = 1.3 MeVd — , (5.29) 



327TaG J Vt ' V t ' V ^ 



so the temperature drops as 1 / v t. 

From this equation one observes that about one microsecond after the Big 
Bang the temperature has dropped from a value above the Planck temperature, 
corresponding to an energy of 10^^ GeV to a temperature of about one GeV, 
so after about one microsecond the temperature is already too low to generate 
protons and after about one second the lightest matter particles, the electrons, 
are "frozen" out. After about three minutes the temperature is so low that the 
light elements become stable, and after about 10^ years atoms can form. 

At this moment all matter becomes neutral and the photons can escape. These 
are the photons, which are still around in the form of the microwave background 
radiation! 

After the discovery of the microwave background radiation, the Big Bang 
theory gained widespread acceptance. Nevertheless, the simplest model as for- 
mulated here, has several serious problems, which can only be solved by the 
so-called inflationary models. These models have the bizarre property that the 
expansion of the universe goes faster than the speed of light! This is not a contra- 
diction of special relativity, since these regions are causally disconnected, so no 
information will be transmitted. Special relativity does not restrict the velocities 
of causally disconnected objects. However, such inflationary scenarios require 
the introduction of a scalar field, e.g. the Higgs field discussed in the previous 
chapter. For certain conditions of the potential of this field, the gravitational 



58 







111)1 
















tot 


\\\\\ 



7 


j 


1 









































































































































































































Figure 5.5: The flatness of tlie universe after inflation is easily understood if one 
tliinks about tlie inflation of a balloon. 



force becomes repulsive, as can be derived directly from the Einstein equations 
given above. This will be discussed in more detail after a short summary of the 
main problems and questions of the simple Big Bang theory. 



5.6 Flatness Problem 

At present we do not know if the universe is open or closed, but experimentally 
the ratio of the actual density to the critical density is bound as follows p[: 

0.1 <n = p/pc< 2. (5.30) 

The luminous matter contributes only about 1% to Q, but from the dynamics of 
the galaxies one estimates that the galaxies contribute between 0.1 and 0.3, so 
the lower limit on Q stems from these observations. The upper limit is obtained 
from the lower limit on the lifetime of the universe. From the dating of the oldest 
stars and the elements one knows that the universe is at least 10^° years old, 
which gives an upper limit on the Hubble constant (eq. |5.19| ) and consequently 
on the density^]. This does look like a perfectly acceptable number and the 
universe might even be perfectly flat, since f2 = 1 is not excluded. However, 
it can be shown easily, that f2 — 1 grows with time as t^^^ and for the present 
lifetime t ^ lO^^s this number becomes big, unless very special initial conditions 
limit the proportionality constant to be exceedingly small. This constant can be 
calculated easily: From eq. p.lO| and the definition of pc (eq. |5.22| ) one finds: 



kc^ 

Since R oc t^^'^ for the longest time of the universe (see table |5.1|) and H(t) oc 1/t 
(eq. |5.21|) one observes that i7 — 1 oc kt'^^^. For this number to come out close 
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Figure 5.6: The space-time diagram of the microwave background radiation, 
which was released at the time tr- Photons observed in point P from opposite 
directions traveled some 10^° yrs with hardly any interactions from the points A 
and B, respectively. The distance light could have traveled between the Big Bang 
and tr is only ab, which is much smaller than the distance AB. Consequently the 
points A and B could never have been in causal contact with each other. Nev- 
ertheless, the radiation from A and B have the same temperature, although the 
horizon ab is much smaller (horizon problem). The problem can be solved if one 
assumes the region of causal contact was much larger than ab through inflation 
of space-time via a phase transition. 



to zero for t very large implies that k must have been very close to zero, right 
from the beginning. Remember that k is proportional to the sum of potential 
and kinetic energy in the non-relativistic approximation. One can showQ] that 
in order for Q to lie in the range close to 1 now, implies that in the early universe 
|fi - 1| < 10-^9 M|,/T2, or for T ^ Mp, 



\n-i\ 



< 10-^^. (5.32) 



In other words, if the density of the initial universe was above the critical density 
say by 10~^^pc, the universe would have collapsed long ago! On the other hand, 
would the density have been below the critical density by a similar amount, the 
present density in the universe would have been negligible small and life could 
not exist! 



5.7 Horizon Problem 

Since the horizon increases linearly with time but the expansion only with t^/^, 
most of the presently visible universe was causally disconnected at the time 
t = 10^ years, when the microwave background was released. Nevertheless, the 
temperature of the microwave background radiation is the same in all directions! 
How did these photons thermalize after being emitted some 10^° years ago? One 
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should realize that the density in the universe is exceedingly low, so photons 
from opposite directions have traveled some 2 ■ 10^° lightyears without interac- 
tions. Since the distance scales as t^/^, these regions were about 10'' light years 
apart at the time they were released, i.e. the distance AB in fig. 5^, which 
is two orders of magnitude larger than the horizon of the universe at that time 
(distance ab), so no signal could have been transmitted. Nevertheless, the tem- 
perature difference AT/T between these regions is less than 10~^ as shown by the 
recent COBE data. As with the flatness problem, one can impose an accidental 
temperature isotropy in the universe as an initial condition, but with all the hefty 
fluctuations during the Big Bang, this is a very unsatisfactory explanation. As 
we will see later, inflation solves both problems in a very elegant way. 



5.8 Magnetic Monopole Problem 

Magnetic monopoles are predicted by GUT's as topological defects in the Higgs 
fleld: after spontaneous symmetry breaking the vacuum obtains a non-zero vac- 
uum expectation value in a given region. Different regions may have different 
orientations of the phases of the Higgs fleld and the borderlines of these regions 
have the properties expected for magnetic monopoles |l55[|. Unfortunately the 



magnetic monopole density is very small, if not zero. Their absence has to be 
explained in any theory based on GUT's with SSB. The flrst attack was made by 
Alan Guth, who invented inflation for this problem. Although the original model 
did not solve the monopole problem, it provided a perfectly reasonable solution 
for the horizon and flatness problem. An alternative version of inflation, the so- 
called new inflation, which was invented by A.D. LindeQ and independently by 
Albrecht and Steinhardt , provided also a solution of the monopole problem. 



as will be discussed in section |5.1CI| . 



5.9 The smoothness Problem 

Our universe has density inhomogeneities in the form of galaxies. On a large scale 
the spectrum of inhomogeneities is approximately scale-invariant, which can be 
understood in the inflationary scenario as follows: the inflation smoothens out 
any inhomogeneities which might have been present in the initial conditions. 
Then in the course of the phase transition inhomogeneities are generated by the 
quantum fluctuations of the Higgs fleld on a very small scale of length, namely 
the scale where quantum effects are important. These density fluctuations are 
then enlarged to an astronomical scale by inflation and they stay scale invariant 
as is obvious if one thinks about a little circle on a balloon, which stays a circle 
after inflation, but just on a larger scale. 



5.10 Inflation 



The deceleration in the universe is given by eq. ^.9| . In case the energy den 



sity only consists of kinetic and gravitational energy, the sign of R is negative. 
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since both the pressure and energy density are positive. However, the situation 
can change drastically, if the universe undergoes a first-order phase transition. In 
Grand Unified Theories such phase transitions are expected: e.g. the highly sym- 
metric phase might have been an SU(5) symmetric state, while the less symmetric 
state corresponds to the SU{3)c ® SU{2)l ® U{1)y symmetry of the Standard 
Model. The description of the spontaneous symmetry breaking by the Higgs 
mechanism leads to a specific picture of this phase-transition: the Higgs field is 
a scalar field, which fills the vacuum with a potential energy V{(j)). The value of 
the potential is temperature dependent: at high temperatures the minimum oc- 
curs for = 0, but for temperatures below the critical temperature, the ground 
state, i.e. the state with the lowest energy, is reached for a value of the field 
7^ 0. This is completely analogous to other phase transitions, e.g. in supercon- 
ductivity the scalar field corresponds to the density of spin Cooper pairs or in 
ferromagnetism it would be the magnetization. 

If during the expansion of the universe the energy density falls below the en- 
ergy density of this scalar field, something dramatic can happen: the deceleration 
can become an acceleration, leading to a rapid expansion of the universe, usually 
called "inflation". 

This can be understood as follows: if the vacuum is filled with this potential 
energy of the scalar field with an energy density pyac, the work W done during 
the expansion is pAV. However, the gain in energy is pyacAV, since the potential 
energy of the vacuum does not change (a "void" stays a "void" as long as no phase 
transition takes place), so an increase in volume implies an increase in energy. 
Since no external energy is supplied, the total energy of the system must stay 
constant, i.e. pAV + p^acAV = 0, or 

P = -pvac- (5.33) 

This is the famous equation of state in case of a potential dominated vacuum. In 
this case eq. |5.1CI| reduces to: 

R = -7iGR{t)p. (5.34) 
3 

This equation has the solution: 

R{t) oc e*/^, (5.35) 

where 

As we have seen in the previous chapter, the symmetry breaking of a GUT 
happens at an energy of 10^^ GeV. The energy density at this energy is extremely 
high: 

u = pc^ = = 10^0° Jm-3, (5.37) 

where the powers are derived from dimensional analysis. Inserting this result into 
eq. |5.36| yields: 

r = lO^^^s. (5.38) 
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Thus the universe inflates extremely rapidly: its diameter doubles every rln2 
s! The behaviour of the Hubble constant during the inflationary era is then 
determined by eq. |5.8|, which yields: 



H(t) = -. (5.39) 

T 

Clearly, the inflationary scenario provides in an elegant way solutions for many 
of the shortcomings of the BBT: 

• The rapid expansion removes all curvature in space-time, thus providing a 
solution for the flatness problem. 

• After expansion the universe reheats because of the quantum fluctuations 
around the new minimum of the vacuum, thus thermalizing a region much 
larger than the visible region at that time. This explains why all visible 
regions in the present universe were in causal contact during the time the 
2.7 K background microwave radiation was released. 

• The rapid expansion explains the absence of magnetic monopoles, since 
after sufficiently large inflation the monopoles are diluted to a negligible 
level. 

However, the whole idea of inflation only works if one assumes the inflation to go 
smoothly from a single homogenous region to a large homogenous region many 
times the size of our universe (so-called new inflation). This requires rather spe- 
cial conditions for the shape of the potential, as was pointed out by Linde0 and 
independently by Albrecht and Steinhardt |^6[ after the original introduction of 
inflation by Guth[Q. The problem is that the corresponding scalar flelds pro- 
viding the potential energy of the vacuum have to be weakly interacting, since 
otherwise the phase transition will involve only microscopic distances according 
to the uncertainty relation. The Higgs flelds providing spontaneous symmetry 
breaking are interacting too strongly, so one has to introduce additional weakly 
interacting scalar flelds. It is non-trivial to combine the requirement of a neg- 
ligible small cosmo logical constant, which represents the potential energy of the 
vacuum, with a vacuum fllled with Higgs fields to generate masses. This re- 
quires large cancellations of positive and negative contributions, which occur e.g. 
in unbroken supersymmetric theories. However, these theories, if they describe 
the real world, have to be broken. Details on these problems are discussed by 
01ive[^ and in the recent text books by Borner|]I[ and Kolb and Turner 0. In 
spite of these problems the arguments in favour of inflation are so strong, that it 
has become the only acceptable paradigm of present cosmology. 

5.11 Origin of Matter 

As discussed above, matter in our universe consists largely of hydrogen (75%) 
and helium (24%) (typically 10^^ nucleons). The absence of antimatter can be 
explained if one has phase transitions and among others CP-violation, as dis- 
cussed in chapter |[ At present the nuclei dominate the energy density in the 
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Figure 5.7: Models including 25% hot dark matter (HDM) can describe the large 
scale structure of the universe, as probed by the galaxy surveys (QDOT) and 
COBE temperature anisotropy, better than models with only dark matter (DM). 
From Schaefer and Shafi [57]. 



universe in contrast to the first 10^ years, when the energy density of the radi- 
ation dominated. The reason for this change is simply the fact that the energy 
density of the many photons around decreases oc T^, while the energy density of 
the nuclei decreases oc (see eq. 

Large amounts of matter can be created from the energy release during the 
inflationary phase of the universe. This can be easily estimated as follows. After 
inflation the universe has a macroscopic size, typically the size of a football or 
larger. The large energy density in this volume (see eq. ^.37| ) yields a total energy 
many times the mass in our present universe. 

Thus within the inflationary scenario the universe could originate as a quan- 
tum fluctuation, starting from absolute "nothing", i.e. a state devoid of space, 
time and matter with a total energy equal to zero. Most matter was created after 
the inflationary phase from the decay of the field quanta of the fields responsable 
for the inflation. Of course, a quantum description of space-time can be discussed 
only in the context of quantum gravity, so these ideas must be considered specu- 
lative until a theory of quantum gravity is formulated and proven by experiment. 
Nevertheless, it is fascinating to contemplate that physical laws may determine 
not only the evolution of our universe, but they may remove also the need for 
assumptions about the initial conditions. 
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5.12 Dark Matter 



The visible matter is clustered in large galaxies, which are themselves clustered in 
clusters and superclusters with immense voids in between. From the movements 
of the galaxies one is forced to conclude that there must be much more matter 
than the observed visible matter, if we want to stick to Newtonian mechanics. 
The most impressive evidence for the dark, i.e. not visible matter comes from the 



so-called flat rotation curves [£8|: the orbital velocities of luminous matter around 
the central of spiral galaxies remain constant out to the far edges of the galaxies 
in apparent contradiction to velocity distributions expected from Keppler's law: 



v\r) = G^^, (5.40) 
r 

where r is the radial distance to the centre of the galaxy, M(r) the mass of the 
galaxy inside a sphere with radius r, and G the gravitational constant. From 
this law one expects the velocities to decrease with if the mass is concen- 

trated in the centre, which is certainly the case for the visible matter. Velocities 
independent of r imply M{r)/r to be constant or M(r) oc r! Such a behaviour 
is expected for weakly interacting matter, like neutrinos, gravitinos or photinos, 
since strongly interacting matter would be attracted to the centre by gravity, 
interact, loose energy, and concentrate in the centre, just like the visible matter 
does. Also the dynamical properties of galaxies in large clusters require large 
amounts of dark matter. To make the speeds work out consistently, one has to 
assume that the total density of 'dark' matter is an order of magnitude more than 
the visible matter. Recent reviews for the experimental evidence of dark matter 



can be found in ref. 59 



From the concentration of light elements, as shown in fig. ^.2| , one has to 
conclude that the total baryonic density is only about 10% of the critical density 
(see eq. p.4| ). The critical density is the density for a flat universe, which naturally 



occurs after inflation. Consequently, in the inflationary scenario the dark matter 
makes up about 90% of the total mass in the universe and it has to be non- 
baryonic. 

Possible candidates for dark matter are the MACHO 's^, which have been 
observed recently through their microlensing effect on the light of stars behind 



them||60|. But since they cluster in heavy compact objects, they are likely to be 
remnants of collapsed stars or light dwarfs, which have too little mass to start 
nuclear burning. In these cases they would be baryonic and make up 10% of the 
critical density, required by baryogenesis. Note that the visible baryonic matter 
in stars represents at most 1% of the critical density. 

So one needs additional dark matter, if one believes in inflation and takes the 
value of fi;, = 0.1 from baryogenesis. Candidates for this additional dark matter 



are neutrinos with a mass in the eV range [pT|: 



9 < < 35 eV. (5.41) 



Such small masses are experimentally not excludedpO], but they would be rela- 



tivistic or "hot". Unfortunately, a// dark matter cannot be relativistic, since this 
''Massive Compact Halo Objects. 
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is inconsistent with the extremely small anisotropy in the microwave background 
as observed by the COBE satellite . 

This anisotropy in the temperature is proportional to the anisotropy in the 
mass density at the time of release of this radiation shortly after the Big Bang. 
Through gravity the galaxies were formed around these fluctuations ( "seeds" ) in 
the mass density. So the present structure of the universe has to follow from 
the spectrum of fluctuations in the early universe, which can be probed by the 
microwave background anisotropy. The best fit is obtained for a mixture of 75% 
cold and 25% hot dark matter|6^, as shown in fig. |5.7| (from ref. iW^l )- 

The lightest supersymmetric particles are ideal candidates for cold dark mat- 
ter, provided they are not too numerous and too heavy. Otherwise they would 
provide a density above the critical density p3|. 



K5| in which case the universe 



would be closed and the lifetime would be very short (see fig. ^.4|) . The LSP can 
annihilate sufficiently rapidly into fermion-antifermion pairs, if the masses of the 
SUSY particles are not too heavy, as will be discussed in chapter 



5.13 Summary 

The Big Bang theory is remarkably successful in explaining the basic observations 
of the universe, i.e. the Hubble expansian, the microwave background radiation 
and the abundance of the elements. From the measured Hubble constant one can 
derive such basic quantities as the the size and the age of the universe. Neverthe- 
less many questions remain unanswered. They can be answered by postulating 
phase transitions during the evolution of the universe from the Planck tempera- 
ture of 10^^ K to the 2.7 K observed today. Among the questions: 

• The Matter- Antimatter Asymmetry in our Universe 

As first spelled out by Sakharov [pT| , any theory trying to explain the pre- 
ponderance of matter in our universe must necessarily implement: 

— Baryon- and Lepton number violation; 

— C- and CP- violation; 

— Thermal non-equilibrium conditions, as expected after phase transi- 
tions. 

• The Dominance of Photons over Baryons 

If the excess of matter originates from small CP-violation effects, most 
matter and antimatter will have enough time to annihilate into photons 
thus providing an explanation why the number of photons as observed in 
the 3K microwave background radiation is about 10^ to 10^° times as high 
as the number of baryons in our universe. 

• Inflation 

An inflationary phase, i.e. a rapid expansion generated by a potential 
energy term, which, according to Einstein's equations of General Relativity, 
provides a repulsive instead of attractive gravitational force, is the only 
viable explanation to solve the following problems: 
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— Horizon Problem 

The fact that the observed temperature of the 2.7 K microwave back- 
ground radiation is to a very high degree the same in all directions can 
only be explained if we assumes that all regions were in causal contact 
with each other at the beginning. However, the size of our universe is 
larger than the "horizon", i.e. the distance light could have traveled 
since the beginning. Therefore, one can only explain the temperature 
isotropy, if one assumes that all regions were in causal contact at the 
beginning and that space expanded faster than the speed of light. This 
is indeed the case in the inflationary scenario. 

— Flatness Problem 

Experimentally the observed density in our universe is close to the 
so-called critical density, which is the density, where the total energy 
of the universe is zero, i.e. the kinetic energy of the expanding uni- 
verse is just compensated by the gravitational potential energy. This 
corresponds to a "fiat" universe, i.e. zero curvature. The inflationary 
scenario naturally explains, why the universe is so fiat: the rapid ex- 
pansion by more than 50 orders of magnitude drives all curvature to 
zeroQ 

— Magnetic Monopole Problem 

Magnetic monopoles are predicted by GUT's. Their absence in our 
universe is explained by the inflationary models, if one assumes the 
inflation to go smoothly from a single homogeneous region to a large 
homogeneous region many times the size of our universe (so-called new 
inflation) . 

— The Smoothness Problem 

Experimentally the cosmic background radiation shows the features in 
accord with the Harrison-Zel'dovich scale-invariant spectrum (n=l)[^, 
which is the spectrum expected after inflation. The scale invariance 
can be understood as follows: the inflation smoothens out out any in- 
homogeneities which might have been present in the initial conditions. 
Then in the course of the phase transition inhomogeneities are gener- 
ated by the quantum fluctuations of the Higgs field on a very small 
scale of length, namely the scale where quantum effects are impor- 
tant. These density fluctuations are then enlarged to an astronomical 
scale by inflation and they stay scale invariant as is obvious if one 
thinks about e.g. a little circle on a balloon, which stays a circle after 
inflation, but just on a larger scale. 

In Grand Unifled Theories the conditions needed for the inflationary Big Bang 
Theory are naturally met: at least two phase transitions, which generate mass 
and thus provide non-equilibrium conditions, are expected: one at the uniflcation 
scale of 10^^ GeV, i.e. at temperature of about 10^^ K and one at the electroweak 
scale, i.e. a temperature of about 10^^ K. Furthermore, a potential energy term in 
the vacuum is expected from the scalar flelds in these theories, which are needed 

^Just like blowing up a balloon removes all wrinkles and curvature from the surface. 
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to generate particle masses in a gauge-invariant way. In the minimal model at 
least 29 scalar fields are required. Unfortunately, none have been discovered so 
far, so little is known about the scalar sector. 

Nevertheless, the arguments in favour of inflation are so strong, that it has 
become the only acceptable paradigm of present cosmology. Experimental obser- 
vation of scalar fields would provide a great boost in the acceptance of the role 
of scalar fields, both in cosmology and particle physics. 
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Chapter 6 



Comparison of GUT's with 
Experimental Data 



In this chapter the various low energy GUT predictions are compared with data. 
The most restrictive constraints are the couphng constant unification combined 
with the lower limits on the proton lifetime. They exclude the SM[^, |n|, 
as well as many other models [O, |35|, B3] with either a more complicated Higgs 



sector or models, in which one searches for the minimum number of new particles 
required to fulfil the constraints mentioned above. From the many models tried, 
only a few yielded unification at the required energies, but these models have 
particles introduced ad-hoc without the appealing properties of Supersymmetry. 
Therefore we will concentrate here on the supersymmetric models and ask if the 
predictions of the simplest, i.e. minimal models|^U[ are consistent with all the 
constraints described in the previous chapters. The relevant RG equations for 
the running of the couplings and the masses are given in the appendix. Assuming 
soft symmetry breaking at the GUT scale, all SUSY masses can be expressed in 
terms of 5 parameters and the masses at low energies are then determined by the 
well known Renormalization Group (RG) equations. So many parameters cannot 
be derived from the unification condition alone. However, further constraints can 
be considered: 



Mz predicted from electroweak symmetry breaking pi, [70| , [fl| , 6^, 72 



b-quark mass predicted from the unification of Yukawa couplings [^, 0, [75 
Constraints from the lower limit on the proton lifetime [ pi] , [7^ , [T^ . 



Constraints on the relic density in the universe [76, 54 



Constraints on the top mass |78|, |70|, |7T], [75 



Experimental lower limits on SUSY masses [^, [79|] . 



Constraints from 6 — > 37 decays[Bn, SI, 82 



Of course, in many of the references given above, several constraints are stud- 
ied simultaneously, since considering one constraint at a time yields only one 
relation between parameters. Trying to find complete solutions with only a few 
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constraints requires then additional assumptions, like naturalness, no-scale mod- 
els, fixed ratios for gaugino- and scalar masses or a fixed ratio for the Higgs 
mixing parameter and the scalar mass, or combinations of these assumptions. 

Several ways to study the constraints simultaneously have been pursued. One 
can either sample the whole parameter space in a systematic or random way and 
check the regions which are allowed by the experimental constraints. 

Alternatively, one can try a statistical analysis, in which all the constraints 
are implemented in a definition and try to find the most probable region of 
the parameter space by minimizing the function. 

In the first case one has to ask: which weight should one give to the various 
regions of parameter space and how large is the parameter space? Some sample 
the space only logarithmically, thus emphasizing the low energy regions[0, oth- 
ers provide a linear sampling [pSj |8^ , |85| , In the second case one is faced with 
the difficulty, that the function to be minimized is not monotonous, because of 
the experimental limits on the particle masses, proton lifetime, relic density and 
so on. At the transitions where these constraints become effective, the derivative 
of the function is not defined. Fortunately, good minimizing programs in mul- 
tidimensional parameter space, which do not rely on the derivatives, exist |]87[ . 
The advantage of such a statistical analysis is that one obtains probabilities for 
the allowed regions of the parameter space and can calculate confidence levels. 
The results of such an analysis p8| will be presented after a short description of 
the experimental input values. Other analysis have obtained similar mass spectra 
for the predicted particles in the MSSM|7y, ^ ^ Q or extended versions 
of the MSSMra. 



6.1 Unification of the Couplings 

In the SM based on the group SU(3) x SU(2) x U(l) the couplings are defined 
as: 

ai = (5/3)^'V(47r) = 5a/(3 cos^ ^vk) 

as = g'^/W) = a/sm^9w (6.1) 

"3 = 9s 

where g' ,g and gs are the f/(l), SU{2) and SU{?>) coupling constants; the first 
two coupling constants are related to the fine structure constant by (see fig. ^.51 ): 



e = v47ra = g sin 6w = g cos Q-w (6.2) 
The factor of 5/3 in the definition of ai has been included for the proper nor- 



malization at the unification point (see eq. p.l4|) . The couplings, when defined 
as effective values including loop corrections in the gauge boson propagators, 
become energy dependent ("running"). A running coupling requires the specifi- 
cation of a renormalization prescription, for which one usually uses the modified 



minimal subtraction (MS*) scheme p7[ . 

In this scheme the world averaged values of the couplings at the Z° energy are 

a-i(Mz) = 127.9 ±0.1 (6.3) 
sin^^^ = 0.2324 ±0.0005 (6.4) 
as = 0.123 ±0.006. (6.5) 
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The value of a ^ is given in ref. ||92| and the value of sin^ Ojja has been been 



taken from a detailed analysis of all available data by Langacker and Polonsky ||93|| , 
which agrees with the latest analysis of the LEP data|2^. The error includes the 



uncertainty from the top quark. We have not used the smaller error of 0.003 for 
a given value of m^, since the fit was only done within the SM, not the MSSM, so 
we prefer to use the more conservative error including the uncertainty from m^. 
The tta value corresponds to the value at Mz as determined from quantities 



calculated in the "Next to Leading Log Approximation" ||94]|. These quantities are 
less sensitive to the renormalization scale, which is an indicator of the unknown 
higher order corrections; they are the dominant uncertainties in quantities relying 
on second order QCD calculations. This Us value is in excellent agreement with 
a preliminary value of 0.120 ± 0.006 from a fit to the cross sections and 



asymmetries measured at LEP|2^, for which the third order QCD corrections 
have been calculated too; the renormalization scale uncertainty is correspondingly 
small. 

The top quark mass was simultaneously fitted and found to be [P^] : 

Mtop = 166^}^ GeV, (6.6) 

where the first error is statistical and the second error corresponds to a variation 
of the Higgs mass between 60 and 1000 GeV. The central value corresponds to a 
Higgs mass of 300 GeV. 

For SUSY models, the dimensional reduction DR scheme is a more appro- 
priate renormalization scheme |^5[. This scheme also has the advantage that all 



thresholds can be treated by simple step approximations. Thus unification occurs 
in the DR scheme if all three meet exactly at one point. This crossing 

point then gives the mass of the heavy gauge bosons. The MS and DR couplings 
differ by a small offset 

^ = (6.7) 

^DR ^MS 127r ^ ' 

where the Ci are the quadratic Casimir coefficients of the group (Cj = N for 
SU(iV) and for U(l) so ai stays the same). Throughout the following, we use 
the 'DR scheme for the MSSM. 



6.2 Mz Constraint from Electroweak Symme- 
try Breaking 

As discussed in chapter ^ the electroweak breaking in the MSSM is triggered by 
the large negative corrections to the mass of one of the Higgs doublets ||34||. After 



including the one-loop corrections to the Higgs potentially, |72|, the following 
expression for Mz can be found (see appendix): 



Ml 



m\ — m\ tan^ [3 — A| 



m 



327r2 M^ 



tan^ 13-1 



2ml + {^l^l - /"^ cot^ (3) 



/Ki) 



^2 



m 



ti 



(6.8) 

SI 



71 



where mi and m2 are the mass parameters in the Higgs potential, tan (3 is the 
mixing angle between the Higgs doublets and the function / has been defined in 
the appendix. The corrections Az are zero if the top- and stop quark masses are 
identical, i.e. if supersymmetry would be exact. They grow with the difference 
— mt"^, so these corrections become unnaturally large for large values of the 
stop masses, as will be discussed later. In addition to relation |6.8| one finds 
from the minimzation of the potential a relation between tan/3 and (see 
appendix), so requiring electroweak breaking effectively reduces the original 5 
free mass parameters to only 3. 



6.3 Evolution of the Masses 

In the soft breaking term of the Lagrangian mo and mi/2 are the universal 
masses of the gauginos and scalar particles at the GUT scale, respectively and 
/i constrains the masses of the Higgsinos. At lower energies the masses of the 
SUSY particles start to differ from these universal masses due to the radiative 
corrections. E.g. the coloured particles get contributions proportional to 
from gluon loops, while the non-coloured ones get contributions depending on 
the electroweak coupling constants only. The evolution of the masses is given by 
the renormalization group equations |]97|, ffT] , which have been summarized in the 
appendix. Approximate numerical mass formulae for the squarks and sleptons, 
mass mixing between the top quarks, gauginos, and the Higgs mass parameters 
are given in chapter |[The exact formulae can be found in the appendix. 



6.4 Proton Lifetime Constraints 

GUT's predict proton decay and the present lower limits on the proton lifetime 
Tp yield quite strong constraints on the GUT scale and the SUSY parameters. As 
mentioned at the beginning, the direct decay p — > e'^Tr^ via s-channel exchange 
requires the GUT scale to be above 10^^ GeV. This is not fulfilled in the Standard 
Model (SM), but always fulfilled in the Minimal Supersymmetric Standard Model 
(MSSM). Therefore we do not consider this constraint. However, the decays via 
box diagrams with winos and Higgsinos predict much shorter lifetimes, especially 
in the preferred mode p —>■ VK^ . From the present experimental lower limit of 
]^g32 y];|3y] for this decay mode Arnowitt and Nath|^ deduce an upper limit on 
the parameter B, which is proportional to l/Tpi 

5 < 293 ± A2{MHj'iMGUT) GeV'^ (6.10) 

Here M^.^^ is the Higgsino mass, which is expected to be above Mqut] else it 
would induce too rapid proton decay. If M^^ would become much larger than 



Mqut, one would enter the non-perturbative regime. Arnowitt and Nath|^8[ give 
the following acceptable range: 

3 < MhJMgut < 10 (6.11) 
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To obtain a conservative upper limit on B, we allow Mh^ to become an order of 
magnitude heavier than Mqut, so we require 



5 < 977 ± 140 GeV^\ (6.12) 

The uncertainties from the unknown heavy Higgs mass are large compared 
with the contributions from the first and third generation, which contribute 
through the mixing in the CKM matrix. Therefore we only consider the sec- 



ond order generation contribution, which can be written as [44 



B = -2(a2/(a3 sin(2/3)) (m^/m|) 10^ GeV-\ (6.13) 



One observes that the upper limit on B favours small gluino masses rrig, large 
squark masses rrig, and small values of tan/3. To fulfil this constraint requires 

tan/? < 10 (6.14) 

for the whole parameter space and requires a minimal value of the parameter mo 
in case mi/2 is not too large, since rUg ^ 2.7mi/2 and m? ~ mo^ + 7mi/2^ (see 
eq. |4.15| ). The constraint can always be fulfilled for very large values of mi/2- 



However, the finetuning constraint implies rrig < 1000 GeV or mi/2 < 350 
GeV. In this case eq. |6.12| requires mo to be above a few hundred GeV, if mi/2 



becomes of the order of 100 GeV or below, as will be discussed below. 



6.5 Top Mass Constraints 

The top mass can be expressed as: 

= (47r)2 Yt{t) sm\(3), (6.15) 



where the running of the Yukawa coupling as function of t = log{^^) is given 
by|3: 

. ^-("ff) . (6.16) 

One observes that Ytit) becomes independent of lt(0) for large values of Yt{Q), 
implying an upper limit on the top mass|]7T], |7S[. Requiring electroweak symmetry 
breaking implies a minimal value of the top Yukawa coupling, typically Yt{Q) > 
(9(10~^). In this case the term 6Yt{0)F{t) in the denominator of |6.16| is much 
larger than one, since F{t) ^ 290 at the weak scale, where t ^ 66. In this case 
Yt{t) = E{t)/QF{t), so from eq. it follows: 



= ^^^^^ sin^(/?) - (190 GeV)^sin^(/5), (6.17) 

where E and F are functions of the couplings only (see appendix). The physical 
(pole) mass is about 6% larger than the running mass[^ |10U[j : 

Mr'' = mt(l + ^ (200 GeV) sin/3, . (6.18) 

V 3 TT / 
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The electroweak breaking conditions require 7r/4</3<7r/2 (eq. [4.45| ); hence 
the equation above imphes for the MSSM approximately: 



145 < Mf'" < 200 GeV, (6.19) 

which is consistent with the experimental value of 166 GeV as determined at 
LEP (see eq. p. 61). Although the latter value was determined from a fit using the 



SM, one does not expect shifts outside the errors, if the fit would be made for 
the MSSM. 

As will be shown in chapter ^ for such large top masses, the b-quark mass be- 
comes a sensitive function of rrit and of the starting values of the gauge couplings 
at Mgut- 

6.6 b-quark Mass Constraint 

As discussed in chapter ^ the masses of the up-type quarks are arbitrary in the 
SU (5) model, but the masses of the down-type quarks are related to the lepton 
masses within a generation, if one assumes unification of the Yukawa couplings at 
the GUT scale. This does not work for the light quarks, but the ratio of b-quark 
and r-lepton masses can be correctly predicted by the radiative corrections to 
the masses [03], [73 |. 



To calculate the experimentally observed mass ratio the second order renor- 
malization group equations for the running masses have to be used. These equa- 
tions are integrated between the value of the physical mass and Mqut- 

For the running mass of the b-quark we used||99||: 



mb = 4.25 ± 0.3 GeV. (6.20) 

This mass depends on the choice of scale and the value of ^^(mb). Consequently, 
we have assigned a rather conservative error of 0.3 GeV instead of the proposed 
value of 0.1 GeV||99[. Note that the running mass (in the MS scheme) is related 
to the physical (pole) mass M^"^*^ by[ |100|| : 



mfe = Mf' f 1 - 1^ _ l2A(^f \ ^ 0.825 M^'^ (6.21) 
so mfe = 4.25 corresponds to M^°'^ ^ 5 GeV. We ignore the running of m-r below 



rrib and use for the pole mass: = 1.7771 ± 0.0005 GeVpUTI. 



6.7 Dark Matter Constraint 

As discussed in chapter |^ there is abundant evidence for the existence of non- 
relativistic, neutral, non-baryonic dark matter in our universe. The lightest su- 
persymmetric particle (LSP) is supposedly stable and would be an ideal candidate 
for dark matter. 

The present lifetime of the universe is at least 10^° years, which implies an 



upper limit on the expansion rate (see eq. |5.19|) and correspondingly on the total 
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relic abundance (compare eq. |5.22|) . Assuming Hq > 0.4 one finds that for eacli 
relic particle species 

n^hl < 1. (6.22) 

This bound can only be obeyed, if most of the LSP's annihilated into fermion- 
antifermion pairs, which in turn would annihilate into photons again. As will be 
shown below, the LSP is most likely a gaugino-like neutralino xo- this case 
the annihilation rate XoXo f f depends most sensitively on the mass of the 
lightest (t-channel) exchanged sfermion: fi^/^o m ^ / ml [p^O^ . Consequently, the 
upper limit on the relic density implies an upper limit on the sfermion mass. 
However, as discussed in chapter the neutralinos are mixtures of gauginos and 
higgsinos. The higgsino component also allows s-channel exchange of the Z° 
and Higgs bosons. The size of the Higgsino component depends on the relative 
sizes of the elements in the mixing matrix |4.3CI| , especially on the mixing angle 



tan P and the size of the parameter /i in comparison to Mi f» 0.4mi/2 and M2 ~ 
0.8mi/2. Consequently, the relic density is a complicated function of the SUSY 
parameters, especially if one takes into account the resonances and thresholds 
in the annihilation cross sections |103||, but in general one finds a large region 



in parameter space where the universe is not overdosed] 76 1. In the preferred 



gaugino-like neutralino region the relic density constraint translates into an upper 



bound of about 1000 GeV on mo |§4[, except for large mi/2, where some SUSY 
masses become much larger than 1 TeV and are therefore disfavoured by the 
fine-tuning criterion (see eq. 



6.8 Experimental lower Limits on SUSY Masses 

SUSY particles have not been found so far and from the searches at LEP one 
knows that the lower limit on the charged leptons and charginos is about half 
the Z° mass (45 GeV)||30[ and the Higgs mass has to be above 62 GeV||79|. The 



lower limit on the lightest neutralino is 18.4 GeV||30|, while the sneutrinos have 
to be above 41 GeV[^. These limits require minimal values for the SUSY mass 
parameters. 



There exist also limits on squark and gluino masses from the hadron colliders 
but these limits depend on the assumed decay modes. Furthermore, if one takes 
the limits given above into account, the constraints from the limits of all other 
particles are usually fulfilled, so they do not provide additional reductions of the 
parameter space in case of the minimal SUSY model. 

6.9 Decay h ^ 

Recently CLEO has published an upper bound for this transition 6 57 < 



5.4 ■ 10 ^[I104||. Furthermore a central value of 3.5 ■ 10 ^ and a lower limit of 



1.5-10 can be extracted from the observed process B -ft"*7||104| and assuming 



that the branching ratio for this process is 15% (using lattice calculations) ||105 

In the SM the transition b —>■ s'j can happen through one-loop diagrams with a 
quark (charge 2/3) and a charged gauge boson. SUSY allows for additional loops 
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involving a charged Higgs and the charginos and neutrahnos[^, |81|, |106| , |82 
In SUSY large cancellations occur, since the W — t and — t loops have an 
opposite sign as compared to the ~ ^ loop and all loops are of the same order 
of magnitude. 

Kane et al. p3 find acceptable rates for the b —>■ sj transition in the MSSM for 



a large range of parameter space, even if they include constraints from electroweak 
symmetry breaking and unification of gauge and Yukawa couplings. They do 
not find as strong lower limits on the charged Higgs boson masses as others H. 



Similar conclusions were reached by Borzumati ||107|| , who used the more complete 
calculations including the flavour changing neutral currents |p.06||. It turns out 



that with the present errors the combination of all constraints discussed above 
are more restrictive than the limits on 6 — 57, so we have not included it in the 
analyses discussed below. 



6.10 Fit Strategy 

As mentioned before, given the five parameters in the MSSM plus acuT and 
Mqut, all other SUSY masses, the b-quark mass, and Mz can be calculated by 
performing the complete evolution of the couplings including all thresholds. 

The proton lifetime prefers small values of tan/5 (eq. p. 14 ), while all SUSY 
masses are expected to be below 1 TeV from the fine-tuning argument (see eq. 



Therefore the following strategy was adopted: uiq and mi/2 were varied be- 
tween and 1000 GeV and tan/5 between 1 and 10. The trilinear coupling ^^(0) 
at Mqut was kept mostly at zero, but the large radiative corrections to it were 
taken into account, so at lower energies it is unequal zero. Varying At{0) between 
+3mo and —Sttiq did not change the results significantly, so the following results 
are quoted for At{0). 

The remaining four parameters - acuT, Mqut, and Yt{0) - were fitted 
with the MINUIT program[87| by minimizing the following function: 



X 



E 



i=l 



(Mz- 91.18)^ 



(mfe-4.25)^ 



{D {771117127713)^ 
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{for D>0) 
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^ {for M > Me^p) 
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The first term is the contribution of the difference between the three calculated 
and measured gauge coupling constants at Mz and the following two terms are 
the contributions from the M^-mass and m^-mass constraints. The last three 
terms impose constraints from the proton lifetime limits, from electroweak sym- 
metry breaking, i.e. D = Vh{vi,V2) — Vh{0,0) < (see eq. [4.51| ), and from 
experimental lower limits on the SUSY masses. The top mass, or equivalently, 
the top Yukawa coupling enters sensitively into the calculation of and Mz- 
Instead of the top Yukawa coupling one could have taken the top mass as a 
parameter. However, if the couplings are evolved from Mqut downwards, it is 
more convenient to run also the Yukawa coupling downward, since the RGE of 
the gauge and Yukawa couplings form a set of coupled differential equations in 
second order (see appendix). Once the Yukawa coupling is known at Mqut, the 
top mass can be calculated at any scale. The top mass can be taken as an input 
parameter too using the value from the LEP data. Unfortunately, the value from 
the LEP data (eq. |2.33|) is not yet very precise compared with the range expected 
in the MSSM (eq. |6.19| ), so it does not provide a sensitive constraint. Instead, we 
prefer to fit the Yukawa coupling in order to obtain the most probable top mass 
in the MSSM. As it turns out, the resulting parameter space from the minimiza- 
tion of this includes the space allowed by the dark matter and the b ^ s'y 
constraints, which have been discussed above. 

The following errors were attributed: o", are the experimental errors in the 
coupling constants, as given above, (jfe=0.3 GeV, 0-5=0. 14 GeV, while an and 
aj^j were set to 10 GeV. The values of the latter errors are not critical, since the 
corresponding terms in the numerator are zero in case of a good fit and even for 
the 90% C.L. limits these constraints could be fulfilled and the was determined 
by the other terms, for which one knows the errors. 

For unification in the DR scheme, all three couplings a~^{fi) must cross at 
a single unification point Mqut UTOSf . Thus in these models one can fit the cou- 



plings at Mz by extrapolating from a single starting point at Mqut back to 
Mz for each of the Ofj's and taking into account all light thresholds. The fitting 
program |]87| will then adjust the starting values of the four high energy parame- 
ters {McuT , c^GUT, and ^^(O)) until the five low energy values (three coupling 
constants, Mz and nib) are "hit". The fit is repeated for all values of mo and 
mi/2 between and 1000 GeV and tan/3 between 1 and 10. Alternatively, fits 
were performed in which mi/2 was left free too. 

The light thresholds are taken into account by changing the coefficients of 
the RGE at the value Q = rrii, where the threshold masses rrii are obtained 
from the analytical solutions of the corresponding RGE (see section [4.4|) . These 
solutions depend on the integration range, which was chosen between rrii and 
Mgut- However, since one does not know mj at the beginning, an iterative 
procedure has to be used: one first uses Mz as a lower integration limit, calculates 
mj, and uses this as lower limit in the next iteration. Of course, since the coupling 
constants are running, the latter have to be iterated too, so the values of aiinii) 
have to be used for calculating the mass at the scale mj|j7y, p.09|| . Usually three 
iterations are enough to find a stable solution. 

Following Ellis, Kelley and Nanopoulos[|66|] the possible effects from heavy 
thresholds are set to zero, since proton lifetime forbids Higgs triplet masses to 
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be below Mqut (see eq. |6.11|) . These heavy thresholds have been considered by 



other authors for different assumptions [ p. 1 Oj , p3| , |111 



SUSY particles influence the evolution only through their appearance in the 
loops, so they enter only in higher order. Therefore it is sufficient to consider 
the loop corrections to the masses in first order, in which case simple analytical 
solutions can be found, even if the one-loop correction to the Higgs potential 
from the top Yukawa couphng is taken into account (see appendix). There is one 
exception: the corrections to the bottom and tau mass are compared directly 
with data, which implies that the second order solutions have to be taken for the 
RGE predicting the ratio of the bottom and tau mass. Since this ratio involves 
the top Yukawa coupling Yf, the RGE for Yf has to be considered in second order 
too. These second order corrections are important for the bottom mass, since the 
strong coupling constant becomes large at the small scale of the bottom mass, 
i.e. as(mfe) ^ 0.2. 

In total one has to solve a system of 18 coupled differential equations: 5 
second order ones(for the 3 gauge couplings, Yt and Yb/Yr) and 13 first order 
ones (for the masses and parameters in the Higgs sector, see appendix). The 
second order ones are solved numerically^ taking into account the thresholds of 
the hght particles using the iteration procedure discussed above. Note that from 
the starting values of all parameters at Mqut one can calculate all light thresholds 
from the simple first order equations before one starts the numerical integration 
of the five second order equations. Consequently, the program is fast in finding 
the optimum solution, even if before each iteration the light thresholds have to 
be recalculated. 



6.11 Results 



The upper part of fig. ^TT] shows the evolution of the coupling constants in 
the MSSM for two cases: one for the minimum value of the function given 



in eq. 6.23| (solid lines) and one corresponding to the 90% CL. upper limit 



of the thresholds of the light SUSY particles (dashed lines). The position of 
the light thresholds is shown in the bottom part as jumps in the first order (3 
coefficients, which are increased according to the entries in table |A.l 



as soon 



as a new threshold is passed. Also the second order coefficients are changed 
correspondingly (see table [A.2| ), but their effect on the evolution is not visible 
in the top figure in contrast to the first order effects, which change the slope of 
the lines considerably in the top figure. One observes that the changes in the 
coupling constants occur in a rather narrow energy regime, so qualitatively this 



picture is very similar to fig. |4.1| , in which case all sparticles were assumed to 
be degenerate at an effective SUSY mass scale M5[/5y[^. Since the running 
of the couplings depends only logarithmically on the sparticle masses, the 90% 
CL. upper limits are as large as several TeV, as shown by the dashed lines in 
fig. ^]l| and more quantitatively in table |6.1| . In this table the initial choices of 



mo and tan/5 as well as the fitted parameters acuT, Mqut, ^^1/2, A*, ^t(O) (and 

^The program DDEQMR from the CERN library was used for the solution of these coupled 
second order differential equations. 
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the corresponding top mass after running down Yt from Mqut to m^) are shown 
at the top and given these parameters the corresponding masses of the SUSY 
particles can be calculated. Their values are given in the lower part of the table. 
Note we fitted here five parameters with five constraints, so the x^=^i if ^ good 
solution can be found. This is indeed the case. The upper and lower limits in 
table |6]l| will be discussed below. 

Only the value of the top Yukawa coupling is given, since for the ratio of 
bottom and tau mass only the ratio of the Yukawa couplings enters, not their ab- 
solute values. For the running of the gauge couplings and the mixing in the quark 
sector, only the small contribution from the top Yukawa coupling is taken into 
account, since for the range of tan/3 considered, all other Yukawa contributions 
are negligible. 

As mentioned before, varying ^^(0) between +3mo and — 3mo does not in- 
fluence the results very much, so its value at the unification scale was kept at 
0, but its non-zero value at lower energies due to the large radiative corrections 
was taken into account. The fits are shown for positive values of the Higgs mix- 
ing paramter /x, but similar values are obtained for negative values of yU with an 
equally good value for the fit. 

The parameters mo, mi/2 and /i are correlated, as shown in fig. |6.2| , where 
the value of /i is shown for all combinations of mg and mi/2 between 100 and 1000 
GeV. One observes that /i increases with increasing itlq and mi/2- The strong 
correlation between mi/2 and /x originates mainly from the electroweak symmetry 
breaking condition, but also from the fact that the thresholds in the running of 
the gauge couplings all have to occur at a similar scale. For example, from fig. 
O it is obvious that the dashed lines for 1/ai and l/a2 will not meet with the 
solid line of l/as, simply because the thresholds are too different; the thresholds 
in 1/0:3 are mainly determined by mi/2, while the thresholds for the upper two 
lines include the winos and higgsinos too, so one obtains automatically a positive 
correlation between and mi/2- 

The value is acceptable in the whole region, except for the regions where 
either itlq or mi/2 or both become very small, as shown in fig. 6^. The increase 
in this corner is completely due to the constraint from the proton lifetime (see 
section |6.5| ). This plot was made for tan/3 = 2. For larger values the region 
excluded by proton decay quickly increases; for tan /? = 10 practically the whole 
region is excluded. 

One notices from fig. |6.2| already a strong correlation between and mi/2- 
This is explicitly shown in fig. |6.4 The steep walls originate from the ex- 
perimental lower limits on the SUSY masses and the requirement of radiative 
symmetry breaking. In the minimum the x^ value is zero, but one notices a long 
valley, where the x^ is only slowly increasing. Consequently, the upper limits on 
the sparticle masses, which grow with increasing values of /i and mi/2, become 
several TeV, as shown in table |6.1| . The 90% CL. upper limits were obtained by 
requiring an increase in x^ of 1.64. so the rather The upper limits are a sensitive 
function of the central value of a^: decreasing the central value of by two 
standard deviations (i.e. 0.012) can increase the thresholds of sparticles several 
TeV. Acceptable fits are only obtained for input as values between 0.108 and 
0.132, if the error is kept at 0.006. Outside this range all requirements cannot be 
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met simultaneously any more, so the MSSM predicts as in this range. 

As discussed previously, sparticle masses in the TeV range spoil the cancella- 
tion of the quadratic divergences. This can be seen explicitly in the corrections 
to Mz- Az is exactly zero if the masses of stop- and top quarks are identical, 
but the corrections grow quickly if the degeneracy is removed, as shown in fig. 
575| . For the SUSY masses at the minimum value of the corrections to Mz are 
small. If one requires that only solutions are allowed for which the corrections 
to Mz are not large compared with Mz itself, one has to limit the mass of the 
heaviest stop quark to about one TeV. The corresponding 90% C.L. upper limits 
of the individual sparticles masses are given in the right hand column of table 
|6.1| . The correction to Mz is 6 times Mz in this case. The limits are obtained 
by scanning tuq and mi/2 till the value increases by 1.64, while optimizing the 
values of tan/3, fi, ,aGUT, ^(0) and Mqut- The lower limits on the SUSY pa- 
rameters are shown in the left column of table |HT^. The lowest values of itlq = 45 
GeV and mi/2=85 GeV are required to have simultaneously a sneutrino mass 
above 42 GeV and a wino mass above 45 GeV. If the proton lifetime is included, 
the minimum value of either mo or mi/2 have to increase (see fig. |6.2|) . Since the 
squarks and gauginos are much more sensitive to mi/2, one obtains the lower lim- 
its by increasing mo. The minimum value for mo is about 400 GeV in this case. 
But in both cases the increase for the lower limits is due to the b-mass, which 
is predicted to be 4.6 GeV from the parameters determining the lower limits. 

The b-quark mass is a strong function of both tan/? and m^, as shown in fig. 
|6.6| ; this dependence originates from the W — t loop to the bottom quark. The 
horizontal band corresponds to the mass of the b-quark after QCD corrections: 
mf, = 4.25 ± 0.3 GeV (see eq. [6.201 ). Since also Mz is a strong function of 
the same parameters, the requirement of gauge and Yukawa coupling unification 
together with electroweak symmetry breaking strongly constrains the SUSY par- 
ticle spectrum. A typical fit with a equal zero is given in the central column 
of table |6.1j , but is should be noted that the values in the other columns provide 
acceptable fits too at the 90% C.L.. 

The mass of the lightest Higgs particle, called h in table 6A., is a rather strong 
function of mt, as shown in fig. for various choices of tan/3, mo and mi/2- All 
other parameters were optimized for these inputs and after the fit the values of 
the Higgs and top mass were calculated and plotted. One observes that the mass 
of the lightest Higgs particle varies between 60 and 150 GeV and the top mass 
between 134 and 190 GeV. Furthermore, it is evident that tan (3 almost uniquely 
determines the value of mt, since even if mi/2 and mo are varied between 100 and 
1000 GeV, one finds the practically the same mt for a given tan/5 and the value 
of mt varies between 134 and 190 GeV, if tan/3 is varied between 1.2 and 5. This 
range is in excellent agreement with the estimates given in eq. |6.19| , if one takes 
into account that Mf"'*^ ^ l.OGmt (see eq. |6.18|). 



Note the strong correlation between tan P and mt in fig. |6.7| : for a given value 
of tan/3 mt is constrained to a vary narrow range almost independent of mi/2 
and mg. Furthermore one observes a rather strong positive correlation between 
mhiggs and all other parameters (tan/?, mo, and mi/2) originating from the 
loop corrections to the potential. 

In summary, the following parameter ranges are allowed (if the limit on proton 
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lifetime is obeyed and extreme finetuning is to be avoided, i.e. mj2 < 1 TeV): 
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The corresponding constraints on the SUSY masses are (see table 3J. for 
details): 
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The lower limits will all increase as soon as the LEP limits on sneutrinos, 
winos and the lightest Higgs increase. 

The lightest Higgs particle is certainly within reach of experiments at present 
or future accelerators [p. 1 2| , 113 |. Its observation in the predicted mass range of 



60 to 150 GeV would be a strong case in support of this minimal version of a 
supersymmetric grand unified theory. 
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Figure 6.1: Evolution of the inverse of the three couphngs in the MSSM. The 
hne above Mqut follows the prediction from the supersymmetric SU(5) model. 
The SUSY thresholds have been indicated in the lower part of the curve: they are 
treated as step functions in the first order j3 coefficients in the renormalization 
group equations, which correspond to a change in slope in the evolution of the 
couplings in the top figure. The dashed lines correspond to the 90% C.L. upper 
limit for the SUSY thresholds. 
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Figure 6.2: The fitted MSSM parameter ji as function of mo and mi/2 for 
tan/3 = 2. 
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Figure 6.3: The of the fit as function of mo and mi/2for tan l3 — 2. The sharp 
increase in in the corner is caused by the lower hmit on the proton hfetime. 
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Figure 6.4: The correlation between mi/2 and jj, for mo=500 GeV. 
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Figure 6.6: The correlation between nib and rrit for mo=400 GeV and two values 
of tan/3. The hatched area indicates the experimental value for m^. 
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Figure 6.7: The mass of the hghtest Higgs particle as function of the top quark 
mass for values of tan/3 between 1.2 and 5 and values of mo and mi/2 between 
100 and 1000 GeV. The parameters of /i, Mqut, oigut and 1^(0) are optimized 
for each choice of these parameters; the corresponding values of the top and 
lightest Higgs mass are shown as symbols. For small values of mi/2 the Higgs 
mass increases with mo, as shown for a "string" of points, each representing a 
step of 100 GeV in mo for a given value of mi/2, which is increasing in steps 
of 100 GeV, starting with the low values for the lowest strings. At high values 
of mi/2 the value of mo becomes irrelevant and the "string" shrinks to a point. 
Note the strong positive correlation between nihiggs and all other parameters: 
the highest value of the Higgs mass corresponds to the maximum values of the 
input parameters, i.e. tan/3 = 5, mo = mi/2 = 1000 GeV; this value does not 
correspond to the minimum y^. More likely values correspond to muggs ~ 92 
GeV for mi/2 = 100, mo = 400 and tan /3 = 2. 
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Tabic 6.1: Values of SUSY masses and parameters for various constraints: 
G=gauge coupling unification; E=electroweak symmetry breaking; Y=Yukawa 
coupling unification; P=Proton lifetime constraint; F=finetuning constraint. 
Constraints in brackets indicate that they are fulfilled but not required. The 
value of the lightest Higgs h can be lower than indicated (see text) . 



Chapter 7 
Summary. 



Many of the questions posed by cosmology suggest phase transitions during the 
evolution of the universe from the Planck temperature of 10^^ K to the 2.7 K 
observed today. Among them the baryon asymmetry in our universe and in- 
flation, which is the only viable solution to explain the horizon problem, the 
flatness problem, the magnetic monopole problem, and the smoothness problem 
(see chapter |^) . 

In Grand Unifled Theories (GUT) phase transitions are expected: one at the 
uniflcation scale of 10^^ GeV, i.e. at a temperature of about 10^^ K and one at 
the electroweak scale, i.e. at a temperature of about 10^"^ K. Furthermore, scalar 
flelds, which are a prerequisite for inflation, are included in GUT's. In the min- 
imal model at least 29 scalar flelds are required. Unfortunately, none have been 
discovered so far, so little is known about the scalar sector, although the verifi- 
cation of the relation between the couplings and the masses of the electroweak 
gauge bosons indeed are indirect evidence that their mass is generated by the 
interaction with a scalar field. Experimental observation of these scalar fields 
would provide a great boost for cosmology and particle physics. First estimates 
of the required mass spectra of the scalar fields can be obtained by comparing 
the experimental consequences of Grand Unified Theories (GUT) with low energy 
phenomenology. 

One of the interesting "discoveries" of LEP was the fact that within the Stan- 
dard Model (SM) unification of the gauge couplings could be excluded (see fig. 
[4.1|) . In contrast, the minimal supersymmetric extension of the SM (MSSM) 
provided perfect unification. This observation boosted the interest in Supersym- 
metry enormously, especially since the MSSM was not "designed" to provide 
unification, but it was invented many years ago and turned out to have very 
interesting properties: 

• Supersymmetry automatically provides gravitational interactions, thus paving 
the road for a "Theory Of Everything" . 

• The symmetry between bosons and fermions alleviates the divergences in 
the radiative corrections, in which case these corrections can be made re- 
sponsible for the electroweak symmetry breaking at a much lower scale than 
the GUT scale. 
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• The lightest supersymmetric partner (LSP) is a natural candidate for non- 
relativistic dark matter in our universe. 

Other non-supersymmetric models can yield unification too,but they do not ex- 
hibit the elegant symmetry properties of supersymmetry, they offer no explana- 
tion for dark matter and no explanation for the electroweak symmetry breaking. 
Furthermore the quadratic divergences in the radiative corrections do not cancel. 

The Minimal Supersymmetric Standard Model (MSSM) model has many pre- 
dictions, which can be compared with experiment, even in the energy range where 
the predicted SUSY particles are out of reach. Among these predictions: 

• Mz. 



• nib. 

• Proton decay. 

• Dark matter. 

It is surprising, that in addition to the unification of the coupling c onstants the 
minimal supersymmetric model can fulfil all experimental constraints from these 
predictions. As far as we know, supersymmetric models are the only ones, which 
are consistent with all these observations simultaneously. Within the MSSM the 
evolution of the universe can be traced back to about 10~^^ seconds after the 



'bang', as sketched in fig. |7]1|. If we believe in the infiationary scenario even 
the actual creation of the universe is describable by physical laws. In this view 
the universe would originate as a quantum fiuctuation, starting from absolute 
"nothing", i.e. a state devoid of space, time and matter with a total energy equal 
to zero. Indeed, estimates of the total positive non-gravitational energy and 
negative potential energy are about equal in our universe, i.e. according to this 
view the universe is the ultimate "free lunch" . All this mass was generated from 
the potential energy of the vacuum, which also caused the infiationary phase. 

Of course, a quantum description of space-time can be discussed only in the 
context of quantum gravity, so these ideas must be considered speculative until 
a renormalizable theory of quantum gravity is formulated and proven by ex- 
periment. Nevertheless, it is fascinating to contemplate that physical laws may 
determine not only the evolution of our universe, but they may remove also the 
need for assumptions about the initial conditions. 

From the experimental constraints at low energies the mass spectra for the 



SUSY particles can be predicted (see table jG^^ in the previous chapter). The 
lightest Higgs particle is certainly within reach of experiments at present or future 
accelerators. Its observation in the predicted mass range of 60 to 150 GeV would 
be a strong case in support of this minimal version of the supersymmetric grand 
unified theory. Discovering also the heavier SUSY particles implies that the 
known strong, electromagnetic and weak forces were all unified into a single 
"primeval" force during the birth of our universe. Future experiments will tell! 
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Figure 7.1: Possible evolution of the radius of the universe and the coupling 
constants. Before t = 10"^^ s spontaneous symmetry breaking occurs, which 
breaks the symmetry of the GUT into the well known symmetries at low energies. 
In the mean time the universe inflates to a size far above the distance light could 
have traveled as indicated by the dashed hne. Prom [114]. 
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Appendix A 



A.l Introduction 

In this appendix all the Renormalization Group Equations (RGE) for the evo- 
lution of the masses and the couplings are given. SUSY particles influence the 
evolution only through their appearance in the loops, so they enter only in higher 
order. Therefore it is sufficient to consider the loop corrections to the masses only 
in first order, in which case a simple analytical solution can be found, even if the 
one-loop correction to the Higgs potential from the top Yukawa coupling is taken 
into account. There is one exception: the corrections to the bottom and tau mass 
are compared directly with data, which implies that the second order solutions 
have to be taken for the RGE predicting the ratio of the bottom and tau mass. 
Since this ratio involves the top Yukawa couphng Y^, the RGE for Yt has to be 
considered in second order too. These second order corrections are important for 
the bottom mass, since the strong coupling constant becomes large at the small 
scale of the bottom mass, i.e. as{mb) ~ 0.2. 

So in total one has to solve a system of 18 coupled differential equations (5 
second order, 13 first order): 

• 3 second order equations for the running of the gauge coupling constants 
ai, i = 1, 3; 

• 2 second order equations for the running of the top Yukawa coupling and 
the ratio of bottom and tau Yukawa coupling Rbr] 

• 1 first order equation for the masses of the left-handed doublet of an u-type 
and d-type squark pair Q; 

• 1 first order equation for the masses of the right-handed up-type squarks 

• 1 first order equation for the masses of the right-handed down-type squarks 
D; 

• 1 first order equation for the masses of the left-handed doublet of sleptons 

• 1 first order equation for the masses of the right-handed singlet of a charged 
lepton E; 

• 4 first order equations for the 4 mass parameters of the Higgs potential 
(mi,m2,m3, and //); 

• 3 first order equations for the (Majorana) masses of the gauginos (Mi,M2 
and Mg). 

• 1 first order equation for the trilinear coupling between left- and right 
handed squarks and the Higgs field At, where the subscript indicates that 
one only considers this coupling for the third generation. 
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Note that the absolute values of the bottom and tau Yukawa couplings need not 
to be known, if one neglects their small contribution to the running of the gauge 
couplings. If one wants to include these, one has to integrate the RGE for Yf, and 
Yt- separately (they are given below too) instead of the RGE for their ratio only. 
Integrating the ratio has the advantage, that the boundary condition at Mqut is 
known to be one, if one assumes Yukawa coupling unification. 

The particle masses are related directly to the Yukawa couplings: 

Ytimt) = j^:^', mt = ht{mt) V sin P (A.l) 

Yb{mb) = j^^; TTib = h{mh) V cos (3 (A. 2) 

F^(mT-) = • rrir = hr{mr) v cos/5. (A.3) 



(47r)2' 



It follows that 



Yi ml 



Yb ml tariff]' ^^'"^^ 

For tan /? < 10, as required by proton decay limits, one observes that Yb is at least 
an order of magnitude smaller than Yf for the values of considered. Hence 
its contribution is indeed negligible in the running of the gauge couplings, so 
below we will only consider the contribution of Yb in the ratio of Yb/Yr, which is 
independent of the absolute value of Yft(O). 

Below we collect all the RGB's in a coherent notation and consider the coef- 
ficients for the various threshold regions, i.e. virtual particles with mass m^ are 
considered to contribute to the running of the gauge coupling constants effec- 
tively only for Q values above mj. Thus the thresholds are treated as simple step 
functions in the coefficients of the RGE. 

Furthermore, all first order solutions are given in an analytical form, including 



the corrections from the top Yukawa coupling ||115|| . Note that a more demanding 
analysis requires a numerical solution of the first five second order equations, in 
which the coefficients are changed according to the thresholds found as analytical 
solutions of the first order equations for the evolution of the masses. The results 
given in chapter |^ all use the numerical solution of these second order equation^. 

Using the supergravity inspired breaking terms, which assume a common mass 
mi/2 for the gauginos and another common mass mo for the scalars, leads to the 
following breaking term in the Lagrangian: 

^Breaking = " "^Q X! I V'* I ^ ~ "^1/2 X! "^Cf-^" (^-5) 

i a 

- Amo [K,QMSH2 + htbQaD'bHi + K^KEIH^] - Bm, [fiH.mQ) 

Here 



^The program DDEQMR from the CERN library was used for the solution of these coupled 
second order differential equations. 
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h^l are the Yukawa couplings, a^h — 1,2,3 run over the generations 

Qa are the SU(2) doublet quark fields 

C/^ are the SU(2) singlet charge-conjugated up-quark fields 

Dl are the SU(2) singlet charge-conjugated down-quark fields 

La are the SU(2) doublet lepton fields 

E'^^ are the SU(2) singlet charge-conjugated lepton fields 

i^i_2 are the SU(2) doublet Higgs fields 

t^i are all scalar fields 

are the gaugino fields 
The last two terms in C Breaking originate from the cubic and quadratic terms 
in the superpotential with A, B and // as free parameters. In total we now have 
three couplings ccj and five mass parameters: 

mo, mi/2, A{t), B{t). 

with the following boundary conditions at Mqut = 0): 

scalars : rfig = = m|, = m\ = m% = ml; (A. 7) 

gauginos : Mi = mi/2, « = 1, 2, 3; (A.8) 
couplings : Oii{0) — (Xgut, ^ = 1,2, 3. (A. 9) 

Here Mi, M2, and M3 are the gauginos masses of the U{1), SU{2) and SU{2>) 
groups. In = 1 supergravity one expects at the Planck scale B = A — 1. 

With these parameters and the initial conditions at the GUT scale the masses 
of all SUSY particles can be calculated via the renormalization group equations. 



A. 2 Gauge Couplings 

The following definitions are used: 

<5i = ^ (A. 10) 

t = In(^) (A.11) 
Pi = hacuT (A.12) 

m = (A.14) 

where 0;^ (i=l,3) denote the three gauge coupling constants of U (1), SU {2)andSU (3), 
respectively, aauT is the common gauge coupling at the GUT scale Mqut and &j 
are the coefficients of the RGE, as defined below. 
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The second order RGE's for the gauge couphngs including the effect of the 



Yukawa couphngs are[0, |7l|, |116 



dai 



where ai = ^, 02 = 6, 03 = 4 for SUSY and ai = j^, a2 
The first order coefficients for the SM are [|117|| : 
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2 for the SM. 



(A.16) 



while for the supersymmetric extension of the SM (to be called MSSM in the 
following) |TT7|: 



(A.17) 



Here Npam is the number of families of matter supermultiplets and Nniggs is 
the number of Higgs doublets. We use Npam = 3 and Nmggs = 1 or 2, which 
corresponds to the minimal SM or minimal SUSY model, respectively. 
The second order coefficients are: 
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For the SUSY model they become: 
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(A.18) 
\ 
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(A.19) 

The contributions for the individual thresholds to hi and hij are listed in tables 
A.l| (from ref. [^) and [A.2| , respectively. 

The running of each ai depends on the values of the two other coupling 
constants, if the second order effects are taken into account. However, these 
effects are small, because the fejj's are multiplied by aj/Air < 0.01. Higher orders 
are presumably even smaller. 

If the small Yukawa couplings are neglected, the RGE's [A.15| can be solved 
by integration to obtain a^(yu') at a scale /i' for a given aj(/i): 



«-(/i') 



/3o-ln 



/2 



2 ai(/i) 




-1 



(A.20) 
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with 



This exact solution to the second order renormahzation group equation can 
be used to calculate the coupling constants at an arbitrary energy, if they have 
been measured at a given energy, i.e. one calculates from a given ai{fi). 

This transcendental equation is most easily solved numerically by iteration. If 
the Yukawa couplings are included, their running has to be considered too and 
one can solve the coupled equations of gauge couplings and Yukawa couplings 
only numerically. 

A. 3 Yukawa Couplings 

In order to calculate the evolution of the Yukawa coupling for the b quark in 
the region between rUb and Mqut, one has to consider four different threshold 
regions: 

• Region I between the typical sparticle masses Msusy and the GUT scale. 

• Region II between Msusy and the top mass m^. 

• Region III between rrit and Mz- 

• Region IV between Mz and m^. 

A. 3.1 RGE for Yukawa Couplings in Region I 

The second order RGE for the three Yukawa couplings of the third generation in 



the regions between Msusy and Mqut are ||118 
dYt ../16^ 13 



Yt —as + 3^2 + —ai - 6Yt 
at V 3 15 

,16, 128,^2 /o, ^^~2 169_2 o~ ~ 136. . 

-16asYt - Q&2Yt - ^a^Yt + 22F/) (A.23) 
dYb ^ /16. ^ 7 . V 



dt V 3 15 

,16,^ 128, ,^2 /o, 9n~2 49_2 o~ ~ 8~ ~ 

~a,Yt + 5yA (A.24) 



dYr , / 9 

Yr +3a2 + -ai 



5 

dt — ■ 5 

2^"' ^5"^ ■ 50'"^ 5 



9 9 81 9 

(362 + - (-61 + —)al - -aatti ) (A.25) 
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If one assumes Yukawa coupling unification for particles belonging to the same 
multiplet, i.e. Yh = at the GUT scale, one can calculate easily the RGE for 
the ratio Rbrif) = mij/mr = ^jYb{t)/YT-{t): 

dRbr 0/8. 2. 1 



dt V3 3 2' 

,8^ 64, 2 /2, 34_2 ~ ^~ ~ 2„ _ 

-(gOs + yjag + + — - 4^302 - -a^ai + -asOi 

hiYt + IyA (A.26) 



5 2 

A. 3. 2 RGE for Yukawa Couplings in Region II 



For the region between Msusy and rrit one finds ||1 19 

dYt ,,(^~ 9. 17. 9,^ 

— — = Yt sa^ H — ao H ai Yt 

dt V 4 20 2 

3 4 2 gQQ 1 3 2 3 1 20 ^ ^ 

22^ 393 \ 

-36(53^; - ^«2l^* - ^aiY, + 12Y,'j (A.27) 



Yft 8^3 + -a2 + -ai - -Y-t 



dt "V 4 ' 4 ' 2' 

2 23 2 127 2 31 27 
+108ao H cto H a, — 90:30:2 03O1 H 02O1 

3 42 1 3 2 15 ^ 20 

99 91 1 \ 

-403^, - ^02^, - ^-a^Y, + -r/ j (A.28) 
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+803^4 &2Yt + OiFt Y; (A.30) 

^ 32 160 4 * / ^ ' 

A. 3. 3 RGE for Yukawa Couplings in Region III 

For the region between mt and Mz one finds: 

dYb ^ / 9 1 
^ = (^803 + -«2 + 

23. 127. . . 31 . . 27 . \ 
+ IO863 + — 02 + - 963O2 - — 0361 + ^"2ai j (A.31) 

dY, /9. 9. 

Yt -Oi2 + -Oi 



dt ' 4 
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23^2 1371^2 27 . \ ,, 

H ttn a, aoQi (A. 32) 

4 2 200 ^ 20 7 ^ ^ 

= Rbr (4^3 - ai 

..~2 53 ^2 9. . 31. . 27. \ 
+54q!3 + — - -0302 - ^"3"! + ^"2ai ) (A.33) 

A. 3. 4 RGE for Yukawa Couplings in Region IV 

dRbr ( ,~ ~ r A ~i 53 _2 31__\ /»oj\ 

= Rbr ( 4^3 - ai + 54^3 + — - ^"3"i 1 (A.34) 

A. 4 Squark and Slepton Masses 

Using the notation introduced at the beginning, the RGE equations for the 



squarks and sleptons can be written as|97 



^ = (aa^Af^ + 2q.iM?) (a,36) 

dTTl^ 16 1 

= (ya3M| + Sa^M^ + — aiM^) - 6,3Yt{ml + rhfj + ml + Alml - /i^) 

(A.37) 

^ = {^-^a.Ml + Y^a.Ml) - 5,,2Y,{ml + ml + + - (A.38) 

The 5j3 factor ensures that this term is only included for the third generation. 

A. 4.1 Solutions for the squark and slepton masses. 

The solutions for the RGE given above are [[7l| : 

= m2 + m?/2«GC/T(^/2(t) + ^/i(t))-cos(2/3)M|(i-sin2^H/) (A.40) 
ml^ = m2 + m?/2«GC/T(^/2(t) + ^/i(t))+cos(2/3)^M| (A.41) 
I + ml/^acuT - cos(2/3)M| sin^ Ow (A.42) 



ml + mj/.acuT (^/al^) + ^/2(t) + ^/i W) " cos(2/3)M|(-l + ^ sin^ 9^) 

(A.43) 



m^^ = + m2/2aGC7T ( ^sit) + ^^(t) + ^/i W) - cos(2/?)M|(i - i sin^ Ow] 



(A.44) 
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m 



ml + mli^acuT ( i^h{t) + — /i(t) ) - cos(2/5)M|(- sin^ Q^] 



(A.45) 

(A.46) 
(A.47) 



For the third generation the effect of the top Yukawa couphng needs to be taken 



into account, in which case the solution given above are changed to [|115|| : 



(A.48) 
(A.49) 



ml-fi^- ml) - ^^GUT (/2(t) + rnl^2 +(Al^O) 



A non-neghgible Yukawa couphng causes a mixing between the weak interaction 
eigenstates. The mass matrix is[pl]: 



MA mo \Hl\+fi\H^,\ 



-ht{At mo \Hl\+^i\Hl\) 
and the mass eigenstates are: 



(A.52) 



tl,2 



+ "^tfl ± V ("^L ~ "^D^ + 4m|( Amo + /i cot 



(A.53) 



A. 5 Higgs Sector 

A. 5.1 Higgs Scalar Potential 



The MSSM has two Higgs doublets (Q = T3 + Yw /2): 

The tree level potential for the neutral sector can be written as: 
V{Hl H^) = ml\H',\' + mim' - m^(i7?i7° + h.c.) + ^^(|i/,T - 

(A.54) 

with the following boundary conditions at the GUT scale mf = m\ = /x^ + 
ml, m\ = — BfimQ, where the value of fi is the one at the GUT scale. 
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The renormalization group equations for the mass parameters in the Higgs 
potential can be written as||97||: 



IT 



3(a2 + -«i - yt)fJ''^ 

3{a2M^ + l&iM^) + 3(a2 + 7"! - Yt)fi'^ 
5 5 



(A.55) 
(A.56) 



^ = 3{a2M^ + laiM^) + 3{a2 + lai)i/ -3Yt{ml+m'lj + ml + A'fml) 
at 5 5 

(A.57) 

drn^ 3 1 1 

= -i&2 + -&i-Yt)ml + 3fxmoYtAt-3fx{a2M2 + -aiMi) (A.58) 

A. 5. 2 Solutions for the Mass Parameters in the Higgs Po- 
tential 



The solutions for the RGE given above are ||97|| : 

fi\t) = g(t)V'(0) (A.59) 

ml{t) = ml + fi\t)+ml^^aGUT{lf2it) + ^h{t)) (A.60) 

ml{t) = q{tffi\0) + ml/^eit) + A(0)momi/2/(t) + ml{h{t) - k{t)Atm%l) 
ml{t) = g(t)m2(0) +r(t)/i(0)mi/2 + s(t)A(0)mo/i(0) (A.62) 



where 



h(t 
kit 

fit 
D{t 

e{t 

s{t 

r{t 
E{t 
Fit 

H2it 



1 



(l + 6Fi(0)F(t))^/^ 
1) 



2'Dit) 
3YtiO)Fit) 
D\t) 

6YtiO)Hsit) 
D\t) 
1 + 6Ft(0)F(t) 

3 \G^it) + y;(0)G2_(t) iH2it) + QYmH,it)f 



3D^it) 



Hsit) 



mO)Fit) 



Dit) 
f3YtiO)Hsit) 



(1 + p^tf^'^^'-^\i + /?2t)'/'^(i + A^)''/^'"'^^ 

Eit')dt' 



16 13 

acuTi—h^it) + 3h2it) + —hiit)) 

6 15 
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H3{t 

Hj{t 

G2{t 
F2{t 



tE{t) - Fit) 
F{t)H2{t) - Hsit) 

aGUT{-^f3{t) + 6Mt)-^Mt)) 



15' 



J Hl{t')E{t')dt' 



aGUT{Sh2{t) + -h{t)) 
5 

aGUT{-lMt) + f2{t)-lm) 
F2it) - ^-Hlit) 

6F3(t) - F,{t) - 4H2it)H^it) + 2F{t)Hl{t) - 2H,{t) 
aGUT{\h{t) + ^h{t)) 

6 15 



F{t)F2{t) - J E{t')F2{t')dt' 



E{t')H^{t')dt' 



The functions fi and hi have been defined before. The Higgs mass spectrum can 
be obtained from the potential given above by diagonahzing the mass matrix: 



1 d^Vn 



(A.63) 



where (pi is a generic notation for the real or imaginary part of the Higgs field. 
Since the Higgs particles are quantum field oscillations around the minimum, 
eq. |A.63| has to be evaluated at the minimum. The mass terms at tree level 



haven been given in the text. However, as discovered a few years ago, the radia- 
tive corrections to the Higgs mass spectrum are not small and one has to take 
the corrections from a heavy top quark into account. In this case the effective 
potential for the neutral sector can be written as [03]: 



2 '2 

,2|tj0|2 I _2|u-0|2 ™2^trOt70 ^ u„ \ \ 3 +9 Qff0|2 [ffO^^^ 



milH'^l' + m'^\H^\' - ml{H^H^ + h.c.) + 



+ 



327r2 



^2 



-4/1 3, -4/1 



^^2 
t2 



A, Tflx <J 

-) - m;(\n —k 



where rhu are field dependent masses, which are obtained from eqns. |A.53| by 
substituting mf = hf 

The minimum of the potential can be found by requiring: 



dV 

am 



2m\vi - 2m\v2 + ^ \^ (f i - vl)vi 
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871- 



;h'ffi{AtmoV2 + fiVi 



(A.64) 



dV 
d\W\ 



2 o2 9 ~^ 9 / 2 2\ 

„, (•^^ -f2)f2 



27712^2 — 2m^Vi — 



h^Atmo{AtmoV2 + fiVi] 



^2 



m 



ti 



+ [(/«) + /(^l) - 2/K)]/.?t;2} = 0, 



where 



m 



f{m ) = m (In — 2 ~ 1 



mi 



From the minimization conditions given above one obtains: 

4 r 



m2tan j3 



(/ + <7^)(tan2/5-l) 

[/(m^J + fi^fnl) - 2/K)] tan^ /5 + (A^m^ tan^ /3 - /x^) 



167r2 



(A.65) 
(A.66) 

(A.67) 

/W) -/(?^?2)1 1 



m^i - m^2 



2mo 



(m^ + m\) sin 2/5 + 



3/1? sin 213 



+ (y4jmo + jjL tan (3) {AtrriQ + /icot /5) 



{/«) + f{ml) - 2/(m?) 



(A.68) 



Ml 



From the above equations one can derive easily: 
m\ — m\ tan^ (3 — A| 



tan^ (3 



S9' 



327r2 cos2 /? 



/«) + /(m^) + 2m,^ + (A,^m^ - /x^ cot^ /3) 



(A.69) 



(A.70) 



Here all rrii are evaluated at Mz using eqns [A.60[ - |A!6^ . Only the splitting 
in the stop sector has been taken into account, since this splitting depends on 
the large Yukawa coupling for the top quark (see the mixing matrix (eq. |A.52| )). 
More general formulae are given in ref. The Higgs masses corresponding to 
this one loop potential are] 



m 



h,H 



ml + ml + A\, 



327r2 sin^ (3 
'327r2 sin^l3Mf 



firhl) + f{ml) + 2m? + (A^m^ + 



^2 \ 

t2J 



™2 



m 



ti 



^^2 



w 



mji - m?2 



m, 



+ M| + An + A 



^22 



(A.71) 
(A.72) 
(A.73) 
(A. 74) 
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An 
A22 



A12 = 



M| + An + A22)' -4m^M| cos^ 2/3 - 4(An A22 - Aj^] 
± -4(cos2 /?M| + sin^ /JM^) A22 -4(sin2 /3M| + cos^ /3M|) An 
\ -4sin2/?(M| + Ml)Ai2 

'"/i(Atmo + cot /S)'' ^ 



3^7 



ml 



167r2 sin2 /^M^, 



167r2 sin^ 



w 



ln( 



4 



^2 



m 



t2 



+ 



2Atmo{Atmo + cot /3) 



~ 1 
m 



ma' 



+ 
^9' 



Atmo^AtniQ + fi cot P) 



™2 



a 



fi{Atmo + /i cot /?) 



167r2 sin^ pM^ 



,7X1 



74tmo(Amo + n cot /?) 



(A.75) 



(A.76) 



m. 



t2 



mil - ml^ 



(A.77) 
(A. 78) 



where 



2n 1 

h{m ) = ^ ^In— , 



d{'m1, ml) —2 — 
and is the mass of a hght squark. 



mf + mo , 
^ ^ In 



2 
mf 



ml 



mi 

2' 

mi 



A. 6 Charginos and Neutralinos 

The RGE group equations for the gaugino masses of the SU{2>), SU{2) and U{1) 
groups are simple: 

^ = -bmM, (A.79) 
with as boundary condition at Mqut- Mi{t = 0) = mi/2- The solutions are: 

(A.80) 



Mi{t) = ^"77^^1/2 



«i(0) 

Since the gluinos obtain corrections from the strong coupling constant a^, they 
grow heavier than the gauginos of the SU{2) group. There is an additional 
complication to calculate the mass eigenstates, since both Higgsinos and gauginos 
are spin 1/2 particles, so the mass eigenstates are in general mixtures of the weak 
interaction eigenstates. 

The mixing of the Higgsinos and gauginos, whose mass eigenstates are called 
charginos and neutralinos for the charged and neutral fields, can be parametrized 
by the following Lagrangian: 

CGaugino-Higgsino = -^MaA^A,, - ^X^^^^X - (V^M^^V + h.C.) 



105 



where Aq, a = 1, 2, . . . , 8, are the Majorana ghiino fields and 



X 



( B \ 
Hi 



V H9 



2 / 




arethe Majorana neutrahno and Dirac chargino fields, respectively. Here all the 
terms in the Lagrangian were assembled into matrix notation (similarly to the 
mass matrix for the mixing between B and in the SM, eq. pT7|) . The mass 



matrices can be written as [|T4| 



Ml 


—Mz cos /3 sin 9w 
\ Mz sin (3 sin 6\y 





Ms 

Mz cos P cos 9w 
—Mz sin P cos 9w 



—Mz cos P sin sin P sin ^ 

M^ cos P cos 6*1^ — M^ sin P cos 
-fi 

-/X y 

(A.81) 

Ma V2Mw sin 

A/2Miy COS P jJ, 

The last matrix has two chargino eigenstates ■Xi2 with mass eigenvalues 



(A.82) 



Ml^ 



Mi + /i' + 2M^ T V (M| - /i2)2 + cos^ 2/5 + 4M^(M| + /i^ + 2M2/i sin 2/?) 

(A.83) 



The four mass eigenstates of the neutralino mass matrix are denoted by x 
1,2,3,4) with masses M^o < ■ ■ ■ < M^o. The sign of the mass eigenvalue 
corresponds to the CP quantum number of the Majorana neutralino state. 

In the limiting case Mi,M2,fi » Mz one can neglect the off-diagonal ele- 
ments and the mass eigenstates become: 



(A.84) 



with eigenvalues |Mi|, IM2I, and respectively. In other words, the bino 
and neutral wino do not mix with each other nor with the Higgsino eigenstates 
in this limiting case. As we will see in a quantitative analysis, the data indeed 
prefers Mi, M2, /i > Mz, so the LSP is bino- like, which has consequences for dark 
matter searches. 



A. 7 RGE for the Trilinear Couplings in the Soft 
Breaking Terms 

The Lagrangian for the soft breaking terms has two free parameters A and B for 
the trilinear coupling and the mixing between the two Higgs doublets, respec- 
tively. 
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Since the A parameter always occurs in conjunction with a Yukawa couphng, 
we will only consider the trilinear coupling for the third generation, called Af. 
The evolutions of the A and B parameters are given by the following RGE||71||: 



dAt Ms M2 13. Ml' 

-03 h 3a2 1" TT'^i — ~ SriAf (A. 85) 



dt V 3 mo mo 15 mo 

dB / M2 1 MA , , 

3 a2— + -«!— - 3FtAt (A.86) 



(it V mo 5 mo 

The B parameter can be replaced by tan (3 through the minimization conditions 
of the potential. The solution for At{t) is: 

^'^'^ - i + 6r,(o)F(t) + ^r^-i + 6r,(o)F(t)J (^-'^^ 



Particle 


hi 


h2 


63 


9 








2 


h 


3 


1 





10 


2 




3 

5 








w 





4 
3 





q-i 


49 

60 


1 


5 
3 


ii 


1 

60 


1 

2 


1 

6 




4 
15 





1 

6 


h 


2 
5 


2 
3 





H 


1 

10 


1 
6 





t 


17 


1 


2 




30 




3 


Standard Model 


41 
10 


19 
6 


-7 


Minimal SUSY 


33 
5 


1 


-3 



Table A.l: Contributions to the first order coefficients of the RGE for the gauge 
coupling constants. 
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Particles 




9 




^ ^ 



^ AS J 




w 




^ ^ 

f 
^ ^ 




qj 


f 19 3 44 \ 

15 5 15 

1 _I 4 

5 3 

1 11 3 _8 1 
\ 30 2 3 / 


Heavy Higgses 
and Higgsinos 




^ A A ^ 

50 10 ^ 

3 29 n 

10 6 

^ ^ 





Table A. 2: Contributions to the second order coefficients of the RGE for the 
gauge couphngs. 
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